An  Investigation  of 
the  Techniques  of  Ellipsometry, 
Internal  Reflection  Spectroscopy,  and 
Moment  Analysis  to  the  Study 
of  Films  and  Substrates 

Kent  F.  Palmer,  Department  of  Physics,  and 
Michael  Z.  Williams,  Department  Df  Mathematics 
Westminster  College,  Fulton,  Missouri 


January  1981 


Final  Report  for  Period  October  1,  1979  -  September  1,  1980 


Approved  for  public  release;  distribution  unlimited. 


•°i  ’'I  !\  S.  .Vs 


ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 


NOTICES 


When  U.  S.  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other 
than  a  definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  nD 
responsibility  nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is 
not  to  be  regarded  by  implication  or  otherwise,  or  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense  Technical  Information  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be  considered  in  any  sense 
as  an  indorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is  releasable  to  the 
National  Technical  Information  Service  (NT1S).  At  NTIS,  it  will  be  available  to  the  general 
public,  including  foreign  nations. 


APPROVAL  STATEMENT 

This  report  has  been  reviewed  and  approved. 


KENNETH  H.  LENERS 
Project  Manager 
Direciorate  of  Technology 


Approved  for  publication: 
FOR  THE  COMMANDER 


MARION  L.  LASTER 
Director  of  Technology 
Deputy  for  Operations 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  HUMBER  2  GOVT  ACCESSION  NO. 

AEDC-TR-80-51 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

AN  INVESTIGATION  OF  THE  TECHNIQUES  OF 
ELLIPSOMETRY ,  INTERNAL  REFLECTION 
SPECTROSCOPY, AND  MOMENT  ANALYSIS  TO  THE 
STUDY  OF  FILMS  AND  SUBSTRATES 

S  TYPE  OF  REPORT  6  PERIOO  C0VERE0 

Final  Report-October  1, 
1979  to  September  1,  1980 

6  PERFORMING  ORG.  REPORT  NUMBER 

7-  AUTHORfl) 

Kent  F.  Palmer  and  Michael  Z.  Williams, 
Westminster  College,  Fulton,  Missouri  65251 

8  CONTRACT  OR  GRANT  NUMBERfa; 

9  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Arnold  Engineering  Development  Center/  DOT 
Air  Force  Systems  Command 

Arnold  Air  Force  Station,  Tennessee  37389 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  9  WORK  UNIT  NUMBERS 

Program  Element  65807F 

II  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Arnold  Engineering  Development  Center/ DOS 
Air  Force  Systems  Command 

Arnold  Air  Force  Station,  Tennessee  37389 

12  REPORT  DATE 

January  1981 

13  NUMBER  OF  PAGES 

69 

14  MONITORING  AGENCY  NAME  «  AOORESSfil  dtlloront  Iron  ConlroUInt  OUjcoJ 

IS  SECURITY  CLASS,  (ol  thla  report) 

UNCLASSIFIED 

15*  DECL  ASSIFlCATlON'OOWNGnAOlNG 
SCHEDULE 

N/A 

16  DISTRIBUTION  STATEMENT  [ol  Ihla  RaporfJ 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  fn  Block  30.  If  different  trom  Report) 


is  supplementary  notes 

Available  in  Defense  Technical  Information  Center  (DTIC) 


19  KEY  WORDS  f  Continue  on  reverse  aide  It  necessary  and  Identify  by  block  number ) 

ellipsometers  substrates  cryogenics 

internal  reflection  infrared  spectra 
spectroscopy  visible  spectra 

measurements  analysis 

films  moments 

20  ABSTRACT  (Continue  on  reverse  aide  it  necessary  end  identity  by  block  number) 

Ellipsometry  and  internal  reflection  spectroscopy  (IRS)  are 
two  methods  of  spectroscopic  measurement  which  have  some  unique 
advantages  over  other  reflection  techniques  for  the  study  of  films 
and  substrates.  Applications  of  these  methods  are  fully  discussed 
for  work  in  the  infrared  and  visible  spectral  regions.  Moment 
analysis  of  spectral  band  shapes  has  been  shown  theoretically  to 
yield  information  on  the  nature  of  the  intermolecular  interactions 


DD  1  JAM  *73  1  473  EDITION  OF  >  NOV  69  IS  OBSOLETE 


UNCLASSIFIED 


UNCLASSIFIED 


20.  ABSTRACT  (Continued) 

in  condensed  phases.  These  data  are  vital  if  concentrations  of 
molecular  species  are  to  be  found  in  a  film.  The  experimental 
difficulties  of  moment  analysis  are  discussed  along  with  some 
preliminary  results  using  existing  AEDC  data.  Areas  for  future 
research  are  suggested  for  the  spectroscopic  analysis  of  films. 


UNCLASSIFIED 


AEDC-TR-80-51 


PREFACE 

The  research  reported  herein  was  performed  by  Westminster  College,  Fulton,  Missouri, 
under  Subcontract  No.  80-25-PWT  for  ARO,  Inc.  (a  Sverdrup  Corporation  Company), 
operating  contractor  of  Arnold  Engineering  Development  Center,  Air  Force  Systems 
Command,  Arnold  Air  Force  Station,  Tennessee.  The  work  covered  the  period  from 
October  1,  1979,  to  September  i,  1980,  and  was  done  under  ARO  Project  Number 
P32K-13C.  Dr.  Herman  E.  Scott,  Directorate  of  Technology,  was  the  Air  Force  project 
manager.  The  manuscript  was  submitted  for  publication  on  September  19,  1980. 


1 


AEDC-TR-80-51 


CONTENTS 

Page 

1.0  INTRODUCTION .  5 

2.0  ELLI PSOMETRY  . 5 

2.1  Reflection  of  a  Monochromatic  Elliptically  Polarized 

Light  Beam .  6 

2.2  The  Ellipsometric  Parameters  ^  and  A .  9 

2.3  Ellipsometer  Components  . 10 

2.3.1  Light  Sources  . 10 

2.3.2  Polarizers  and  Analyzers . 10 

2.3.3  Compensators  or  Retarders . 11 

2.3.4  Detector  Systems  . 13 

2.4  Ellipsometric  Techniques  for  the  Measurements  of  ^  and  A  . 13 

2.4.1  The  Variable  Retardation  Null  Method . 13 

2.4.2  The  Quarter-Wavelength  Compensator  Null  Method . 14 

2.4.3  The  Fixed  Compensator  Null  Method  . 17 

2.4.4  Principal  Angle  Null  Methods  . 23 

2.4.5  Nonextinction  Methods . 24 

2.5  Errors  Associated  with  Ellipsometer  Components  . 25 

2.6  Relationship  of  the* Ellipsometric  Parameters  ^  and  A 

to  the  Optical  Properties  of  the  Sample  System . 25 

2.6.1  The  Ambient -Substrate  System  . 27 

2.6.2  The  Ambient-Thin  Film-Substrate  System  . 28 

2.7  Commercially  Available  Ellipsometric  Instrumentation  . 31 

2.7.1  Commercial  Production  Ellipsometers  . 32 

2.7.2  Commercial  Research  Ellipsometers  . 33 

2.7.3  Recommendations  Concerning  Ellipsometric  Equipment . 34 

2.8  Summary  . 35 

3.0  INTERNAL  REFLECTION  SPECTROSCOPY  (IRS)  OF  FILMS  . 36 

3.1  Theory  of  Total  Reflection  at  the  Interface  of  Two 

Semi-Infinite  Media  . 37 

3.1.1  Total  Reflection  for  Two  Transparent  Media . 37 

3.1.2  Attenuated  Total  Reflection  for  a  Semi-Infinite 
Transparent  Dense  Medium  and  a  Semi-Infinite 

Absorbing  Medium  . 46 

3.1.3  Effective  Thicknesses  of  a  Semi-Infinite  Weakly 

Absorbing  Rarer  Medium  . 47 

3.1.4  Effective  Thicknesses  of  a  Very  Thin  Weakly 

Absorbing  Film . 48 


3 


A  E  DC-TR  -80-51 


3.2  Experimental  Techniques . 49 

3.1.1  Internal  Reflection  Elements . 50 

3.2.2  Single  Reflection  Methods . 50 

3.2.3  Multiple  Internal  Reflection  (MIR)  Techniques  . 51 

3.3  Commercial  Instrumentation  for  Internal  Reflection 

Spectroscopy  (IRS)  . 53 

3.4  Summary  on  Internal  Reflection  Spectroscopy  of  Films  . 54 

4.0  ANALYSES  OF  MOLECULAR  ABSORPTION  BANDS  BY  THE 

METHOD  OF  MOMENTS  . 54 

4.1  The  First  and  Second  Moments . 54 

4.2  Computer  Program  MOMENT . 55 

4.3  Analyses  with  Program  MOMENT  . 56 

4.4  Mathematical  Band  Models . 57 

4.5  Comparison  with  Other  Moment  Analyses  . 58 

4.6  Summary  on  Moments . 59 

5.0  SUMMARY  . ' . 59 

6.0  REFERENCES . 60 


ILLUSTRATIONS 


Figure 

1 .  Reflection  of  an  Elliptically  Polarized  Light  Wave  .  7 

2.  Typical  Ellipsometer  Arrangement:  the  Polarizer-Compensator-Sample 

System-Analyzer  (PCSA)  Configuration . 12 

3.  Electric  Field  of  a  Beam  Emerging  from  the  Polarizer . 17 

4.  Electric  Field  of  the  Plane-Polarized  Beam  Incident  upon  the  Analyzer . 22 

5.  Relationship  of  the  Two  Possible  Polarizer  Angle  Values,  P[  and  Pj  . 23 

6.  Reflection  and  Refraction  for  an  Ambient-Substrate  System . 26 

7.  Reflection  and  Refraction  for  an  Ambient-Thin  Film- Substrate  System  . 29 

8.  Parallel  and  Perpendicular  Components  of  the  Electric  Field  of  the 
Incident  (E),  Reflected  (R),  and  Refracted  (T)  Beams  at  the  Interface 
of  Two  Transparent  Media  when  the  Incident  Beam  is  in  the  Optically 

Denser  Medium . 39 


TABLE 


1 .  Commercially  Available  Ellipsometers  . 31 

APPENDIX  . 67 


4 


AEDC-TR-80-51 


1.0  INTRODUCTION 

Spectroscopic  techniques  have  been  used  extensively  for  the  identification  of  molecular 
species  and  the  determination  of  molecular  structure.  For  gases  and  condensed  phases  whose 
absorption  laws  are  known,  molecular  species  concentrations  can  be  found.  The  accuracies 
of  molecular  identification  and  species  concentration  depend  on  the  precision  of  the 
spectroscopic  methods.  In  this  report  two  reflection  sepctroscopic  methods  prove,  in  some 
instances,  to  be  superior  to  more  conventional  reflection  or  transmission  techniques. 

At  the  turn  of  the  century,  Paul  Drude  developed  a  method  for  the  analysis  of  polarized 
light  reflected  from  surfaces.  The  method,  now  called  ellipsometry,  is  especially  sensitive  to 
thin-film  formation  on  a  substrate,  but  it  is  primarily  used  for  the  near  infrared  (IR>,  visible, 
and  near  ultraviolet  (UV)  spectral  regions  where  suitable  optics  are  readily  available.  On  the 
other  hand,  internal  reflection  spectroscopy  (IRS)  can  be  readily  applied  throughout  the  IR 
spectral  region.  IRS  spectra  are  very  sensitive  to  the  optical  constants,  and  even  weak 
absorptions  can  be  observed  by  using  multiple  reflections.  Ellipsometry  and  IRS  methods 
are  explored  in  Sections  2.0  and  3.0. 

In  condensed  phases  of  molecules,  intermolecular  effects  become  much  more  important 
than  they  are  for  the  gas  phase,  and  condensed  phase  spectra  have  broader  and  more 
blended  absorption  features  than  do  the  gas  phase  spectra  of  the  same  species. 
Consideration  of  the  molecular  interactions  is  important  if  molecular  species  concentrations 
are  to  be  found  in  solid  or  liquid  films  (Refs.  1  and  2).  Because  moment  analyses  of  spectral 
band  shape  can  be  interpreted  in  terms  of  intermolecular  interactions,  we  have  thoroughly 
investigated  the  band  shapes  of  existing  AEDC  cryofilm  data  and  presented  the  results  in 
Section  4.0.  Areas  for  future  research  are  suggested  in  Section  5.0. 

2.0  ELLIPSOMETRY 

A  light  beam  reflected  from  an  interface  of  two  media  is  altered  in  its  state  of 
polarization  from  the  incident  beam.  The  change  in  polarization  is  affected  by  the  optical 

i 

constants  and  the  geometry  of  the  media.  If,  in  addition,  other  nearby  interfaces  can  cause 
multiple  reflection  effects  (e.g.,  the  two  interfaces  of  a  thin  lamellate  film  sandwiched 
between  a  semi-infinite  ambient  medium  and  a  semi-infinite  substrate),  then  the  optical 
constants  and  geometry  of  all  interacting  media  affect  the  change  in  polarization.  Analyses 
of  the  changes  in  polarization  by  the  techniques  of  ellipsometry  are  especially  useful  for 
determining  the  optical  properties  of  a  substrate  material  having  a  large  extinction 
coefficient,  k,  or  finding  the  optical  constants  and  thickness  of  a  very  thin,  uniform  film  on 
a  substrate.  The  use  of  small-diameter  light  beams  allows  the  extraction  of  spatial 
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information  about  surfaces.  However,  the  standard  optical  materials  used  in  currently 
available  commercial  ellipsometers  limit  the  spectral  domain  to  the  visible,  near-IR,  and 
near-UV.  With  optional  quartz  optics,  measurements  may  be  made  from  2.5  micrometers  to 
0.21  /im.  Ellipsometric  methods  in  the  IR  will  be  discussed  later. 

An  article  by  Neal  (Ref.  3)  contains  an  excellent  overview  of  ellipsometry  and  an 
extensive  bibliography.  A  more  advanced  and  thorough  treatment  of  ellipsometric  methods 
is  given  by  Assam  and  Bashara  (Ref.  4). 

2.1  REFLECTION  OF  A  MONOCHROMATIC 
ELLIPTICALLY  POLARIZED  LIGHT  BEAM 

1 

A  classical  monochromatic  beam  has  a  wavetrain  that  is  infinite  in  extent  and  completely 
polarized  (ReL  5).  The  time-  and  space-variant  electric  field  associated  with  the  light  beam 
can  be  represented  as  the  resultant  of  two  mutually  perpendicular  vectors  having  a  constant 
phase  difference  between  them. 

Let  the  electric  field  disturbance  E  of  a  beam  propagating  in  the  z  direction  of  a  right- 
handed  coordinate  system  be  represented  by  the  real  part  of  E(z,  t)u,  where  the  circumflex 
denotes  a  complex  quantity,  t  is  time,  and  ti  is  a  unit  vector  in  the  x-y  plane.  In  Fig.  1 
(adapted  from  Ref.  6)  x  is  oriented  parallel  to  the  plane  of  incidence.  fi(z,  t)u  can  be 
decomposed  into  the  sum  of  the  p  component 

0  =  |rp«exP  +  0) 

which  is  parallel  to  the  plane  of  incidence  in  the  direction  of  unit  vector  up  and  the  s 
component 

Es  0  =  Lso  op  j-  i  [<a  (L  -  n0  z/c)  +  as  ]  j  (2) 

which  is  perpendicular  to  the  plane  of  incidence  in  the  direction  of  the  unit  vector  us. 
Thus, 

E(*,t)ff  =  Es  U,  0  us  +  Ep  (?,  t)  up 

=  exp  j-  i  (t  -  n0  z/c  )  +  a^j|up,f=p,  s  (3) 

Here,  cj  =  2ircv  is  the  angular  frequency  of  the  wave,  c  is  the  speed  of  light  in  vacuo,  and  v  is 
the  vacuum  wavenumber,  which  is  the  reciprocal  of  the  vacuum  wavelength,  X.  Note  that  in 
Muller’s  convention  for  ellipsometric  notation  (Ref.  7),  the  time  dependence  of  the 
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electromagnetic  field  is  exp(  +  iwt)  so  that  the  complex  index  of  refraction  of  a  medium  is 
n  =  n  -  iK.  In  spectroscopy,  the  time  dependence  is  usually  as  exp(— iwt)  so  that 
n  =  n  +  iK, 

In  Fig.  1,  the  constant  phase  difference  a  =  av  -  as  is  indicated  as  well  as  the  relative 
positions  of  the  maxima  of  each  component.  In  all  practical  cases  the  incident  medium  is 
nonabsorbing  (i.e.,  ri0  =  n0).  The  wavetrain  is  plane-polarized  if  a  is  an  integer  multiple  of 
x,  and  elliptically  polarized  otherwise. 

z" 


Figure  1.  Reflection  of  an  elliptically  polarized 
light  wave. 

From  Fig.  1  and  Eqs.  (1)  through  (3),  we  may  write 

F.p  (z,  t)  =  Re  | Ep  (z,  t)  |  =  Epocos(0  +  a) 

and 

Es  (z,  t)  =  Re  |  E^  (z,  t)  |  =  Eso  cob  6 

for  the  p  and  s  components  of  the  electric  field,  where 

8  =  w  (t  -  n0  z/c )  +  as 

is  a  function sof  t  and  z.  From  Eq.  (4)  we  have 

Ep/Epo-C0s5  cos o  -  -sinS  sina 


(4) 

(5) 

(6) 
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which  upon  squaring  yields 

^E2/E2\-2(E  /E  1  cos  6  cos  a  +  cos2  8  cos2  a  =  sin 2  {?  sin* a 

l  pf  pol  p/  po 

Substituting  cos  6  =  E5/Eso  from  Eq.  (5),  then  rearranging  terms,  gives 


E2  E2  2  E  E  oqs a 

p  B  P  e 


sinZa 


(7) 


Equation  (7)  describes  an  ellipse,  as  can  be  shown  by  the  descriminant  test  described  in 
analytical  geometry  texts.  The  semimajor  axis,  a,  and  semiminor  axis,  b,  are 


These  results  may  be  ascertained  by  performing  a  rotation  by  ij  of  the  x-  and  y-axes  about  the 
z-axis  of  Fig.  1  with  use  of  the  transformation 

•■in 

COS  Tt/  \  E^/ 

_  —  A 

where  E*  and  Ey  are  the  components  of  E  =  Re(E)u  that  are  parallel  to  the  semimajor 
and  semiminor  axes,  respectively. 


The  value  of  the  phase  difference  a  between  the  p  and  s  components  determines  not  only 
the  shape  and  orientation  of  the  elliptical  path  traced  by  the  electric  field  vector  E,  as  the 
above  equations  indicate,  but  also  the  sense  of  its  rotation.  If  a  is  an  integer  multiple  of  x  , 
Eq.  (7)  reduces  to  a  constant  proportionality  between  Ep  and  Es,  i.e.. 


and  the  orientation  of  E  is  confined  to  a  plane,  i.e.,  U  is  constant.  If  a  is  within  the  first 
or  second  quadrants,  E  appears  to  rotate  clockwise  when  looking  into  the  beam.  If  a  is 
within  the  third  or  fourth  quadrants,  the  rotation  is  counterclockwise. 
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Upon  reflection  at  the  interface  of  two  media,  the  electric  field,  R,  of  the  reflected 
beam  can  be  represented  by 


R  (z t)  =  Re 


where 


{r  (*',  t)}  S' 


(ID 


ru;  oS'  =  r  (z',  1)  s '  +  Rfz'.os 

r  F  &  s 


=  SRe0exp|-i[fl»(t-n(|zVc)  t  /S£]j  Sf 


f  =  p,  s 

Here,  z’  is  the  direction  of  propagation  of  the  reflected  wave,  U'  is  in  the  x'-y'  plane  of 
the  x'-  y'-  z'  axes  of  Fig.  1,  U'p  and  U',  (both  perpendicular  to  z')  are  unit  vectors1 
along  the  p  and^ directions,  respectively,  the  constant  phase  difference  between  the  p  and  s 
components  of  R  is  jS  =  j3p  -  /3S,  and 

Rp(z,  1)  =  Rgo  exp|-i[<u  (t  -nQz'/c)  -  jS£]  j 

£  =  P.  s  (12) 

In  general,  R  is  also  elliptically  polarized:  Again,  note  that  n0  is  usually  real  in 
ellipsometric  experiments. 

2.2  THE  ELLIPSOMETRIC  PARAMETERS  ^  and  A 


An  ellipsometer  measures  the  change  in  the  polarization  states  of  the  incident  and  the 
reflected  (or,  very  rarely,  the  refracted)  light  beams.  In  essence,  Es  is  altered  in  the  reflection 
process  to  become  R*.  and  Ep  changes  to  Rp.  The  information  about  the  change  in 
polarization  is  contained  in  the  complex  ratio  of  ratios 


R  /  E  H  /It 

p /  P  _  p/  E 

R  /E  E  /E. 

S,1  F  p;  5- 


(13) 


evaluated  at  the  boundary.  Equations  (1),  (2),  (12),  and  (13)  give 


R  /  R 

po/  so 

exp  | 

0.  -*.>] 

E  / F'. 

po/  so 

exp 

[i(a.-“P>] 

(14) 


=  tan  t b  exp  (-i  A) 
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where 


tnn  i b 


IR  /R 


E  /E 

po  /  so 


(15) 


and 


A  =  -  a  (16) 

are  the  conventionally  defined  quantities  for  the  ratio  of  the  electric  field  amplitude  ratios 
and  the  change  in  the  phase  difference,  respectively. 


2.3  EL1PSOMETER  COMPONENTS 


An  ellipsometer  is  an  instrument  that  can  produce  a  collimated  beam  of  elliptically 
polarized  light  which  can  be  reflected  from  a  sample  system  (S)  and  then  checked  for  the 
changes  of  polarization  caused  by  the  reflection  process.  Many  instrument  geometries  are 
possible,  but  all  ellipsometers  include  (1)  a  collimated  monochromatic  light  source,  L;  (2)  a 
polarizer,  P;  (3)  a  sample  system  mount,  which  may  be  oriented  vertically  (for  solid  samples 
only)  or  horizontally  (for  solid  and  liquid  samples);  (4)  an  analyzer,  A,  itself  a  polarizer,  and 
(5)  a  detector  system,  D.  In  addition,  most  ellipsometers  contain  (6)  a  compensator,  C. 
Figure  2,  adapted  from  Refs.  4  and  6,  shows  a  typical  ellipsometer  geometry  known  as  the 
PCS  A  configuration. 

2.3.1  Light  Sources 

The  collimated  monochromatic  light  source  may  consist  of  a  conventional  polychromatic 
light  source,  a  narrow  band  filter,  slit,  and  collimating  lens  system.  Lasers  are  also  used. 
Commercial  optics  are  available  that  allow  a  He-Ne  (helium-neon)  laser  beam  (0.6328  fim)  to 
be  reduced  to  a  cross-sectional  diameter  of  25  /im  or  smaller.  This  allows  spatial  testing  of 
thin  films  and  surfaces. 


2.3.2  Polarizers  and  Analyzers 

The  partially  polarized  or  unpolarized  radiation  from  a  light  source  will  be  plane- 
polarized  upon  passing  through  a  polarizer.  Nicol,  Gian,  and  Glan-Thompson  prisms  are 
often  used  as  polarizers  or  analyzers  from  the  near  IR  to  the  near  UV.  In  the  medium  and  far 
IR,  the  non-normal  transmission  of  a  beam  through  stacks  of  thin  transparent  plates  can 
cause  satisfactory  plane-polarization.  Other  IR  polarizers  include  the  double  Gian  prism 
(Ref.  8).  Plane-polarized  light  from  the  polarizer  can  be  converted  to  elliptically-polarized 
light  by  a  compensator.  The  orientation  of  a  plane-polarized  wave  can  be  found  with  the 
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orientation  of  a  polarizer  used  as  an  analyzer.  The  polarizer,  analyzer,  and  compensator  are 
mounted  on  rotators  capable  of  accurate  angular  determination,  e.g.,  ±0.01  deg. 

2.3,3  Compensators  or  Retarders 

The  phase  difference  between  two  components  of  a  wave  can  be  altered  upon  passing 
through  a  compensator,  which  allows  one  component  to  pass  through  faster  than  the  other. 
Mica  and  calcite  are  two  examples  of  noncubic  anisotopic  crystalline  materials  which  are 
birefringent,  i.e.,  materials  within  which  a  normally  incident  beam  of  light  is  divided  into 
two  rays,  (1)  a  ray  nearly  undeviated  in  its  path  (the  “ordinary”  ray)  and  (2)  an 
“extraordinary”  ray,  which  is  noticeably  deflected.  A  thin  birefringent  plate  will  allow  a 
normally  incident  plane-polarized  wave  to  exit  a  plane-polarized  wave  only  when  its  plane  of 
polarization  cuts  the  plate  surfaces  in  two  mutually  perpendicular  directions  in  the  plane  of 
the  plate  surfaces. 

The  directions  of  the  lines  of  intersection  are  the  privileged  directions;  the  speeds  of 
propagation  differ  in  the  two  privileged  directions.  In  the  plate,  the  “slow”  component  of 
an  incident  wave  parallel  to  the  slow  direction  follows  the  path  of  the  ordinary  ray  with  an 
effective  index  n0,  and  the  “fast”  component  follows  the  extraordinary  path  with  an 
effective  index  ne  <  no. 

The  phase  shift  6,  or  the  relative  retardation,  of  the  plate  can  be  found  by  considering  the 
fast  component,  E0(z,  t),  and  the  slow  component,  Ee(z,  t),  of  the  electric  vector,  E(z,  t), 
of  a  plane-polarized  light  beam  normally  incident  upon  a  thin  birefringent  plate.  If  z  =  o  at 
the  front  surface,  then  the  form  of  each  component  within  the  plate  is,  from  Eqs.  (1)  and 
(2), 


If  the  deflection  of  the  two  components  is  negligible  for  a  thin  birefringent  plate  of  thickness 
r,  the  final  change  in  phase  between  components  is 

s  =  ~T  (n0-nc)  =  2™7-(no-ne)  (17) 

If  ellipsometric  measurements  are  made  at  a  single  vacuum  wavenumber,  v,  the  equations 
given  below  can  be  simplified  if  r  is  chosen  so  that  5  =  n/2.  Such  a  plate  is  termed  a  quarter- 
wave  compensator. 
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Phase  differences,  6,  can  be  varied  during  an  experiment  by  using  a  Babinet-Soleil 
compensator  (Ref.  5),  or  a  Pockels  effect  retarder  (Ref.  9).  The  latter  consists  of  a  plate  of 
ammonium  dihydrogen  phosphate  (ADP),  which  becomes  birefringent  when  a  potential 
difference  exists  across  its  thickness. 

2.3.4  Detector  Systems 

The  detector  system  is  usually  an  electro-optical  device  such  as  a  photomultiplier  tube 
(PMT)  for  near  IR  to  near  UV  measurements,  or,  for  1R  work,  a  photovoltaic, 
photoconductor,  or  pyroelectric  detector.  For  measurements  over  a  wide  spectral  range,  the 
detector  system  can  be  a  spectrophotometer  (Ref.  8).  Often  a  telescope  is  used  in 
conjunction  with  the  electro-optical  system  to  make  the  preliminary  settings  of  the  optical 
components.  Fine  adjustments  are  made  either  manually  or  automatically  by  a  feedback 
signal  monitoring  the  output  of  the  detector  system. 

2.4  ELLIPSOMETRIC  TECHNIQUES  FOR  THE 
MEASUREMENTS  OF  \p  AND  A 

The  methods  of  eilipsometrie  measurements  fall  into  five  categories  (Ref.  10).  The  first 
four  discussed  are  null  techniques  in  which  ellipsometer  components  are  rotated  to  allow  a 
minimum  of  light  to  pass  through  an  analyzer.  Thus,  the  light  incident  upon  the  analyzer 
must  be  plane-polarized,  or  very  nearly  so. 

If  the  compensator  precedes  the  sample  system  in  the  optical  train,  as  in  the  PCSA 
arrangement  of  Fig.  2,  the  orientations  of  the  components  must  be  such  that  the  beam 
incident  upon  the  sample  is  properly  elliptically  polarized  so  that  the  reflected  beam  is  plane- 
polarized.  (In  the  PSCA  arrangement,  the  compensator  follows  the  sample  system,  the  beam 
incident  upon  the  sample  is  plane-polarized.)  The  elliptically-polarized  beam  reflected  from 
the  sample  becomes  plane-polarized  again  after  passing  through  the  properly  oriented 
compensator.  The  four  general  methods  of  null  ellipsometry  differ  in  regard  to  which 
components  are  adjusted  to  achieve  the  null  condition,  namely  the  extinction  of  the  beam  by 
the  analyzer. 

2.4.1  The  Variable  Retardation  Null  Method 

The  variable  retardation  method  permits  a  direct  measurement  of  the  phase  difference, 
A.  At  a  fixed  orientation  of  the  fast  axis,  the  phase  shift,  6,  of  a  variable  compensator  (e.g., 
a  Babinet-Soleil  compensator  or  a  Pockels  effect  retarder)  is  adjusted  as  well  as  the 
orientation  of  either  the  polarizer  or  analyzer.  At  null,  the  phase  difference  of  the  beam  and 
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the  phase  shift  of  the  compensator  are  equal  (i.e.,  A  =  6),  and  the  orientation  of  the  rotating 
polarizer,  or  analyzer,  gives  ip. 

This  method  is  not  as  accurate  as  other  methods  if  manual  adjustments  of  5  are  made. 
However,  the  use  of  an  electro-optical  retarder  gives  a  measuring  precision  comparable  to 
that  of  other  null  methods. 

2.4.2  The  Quarter-Wavelength  Compensator  Null  Method 

In  the  quarter-wavelength  compensator  method,  either  the  polarizer  (in  the  PSCA 
configuration)  or  the  analyzer  (in  a  PCSA  arrangement)  has  its  axis  set  at  ±  (ir/4)  relative  to 
the  plane  of  incidence.  The  unfixed  component  and  compensator  are  rotated  until  the  null 
condition  appears.  Some  care  must  be  exercised  because  the  settings  of  the  two  components 
are  not  independently  adjustable. 

Deriving  the  appropriate  relations  for  computing  the  ^  and  A  parameters  from  the 
experimental  settings  involves  much  algebra  (Ref.  11).  If  in  a  PSCA  arrangement,  the 
polarizer  axis  setting  is  ±  (tt/4)  from  the  plane  of  incidence,  if  the  angle  of  the  fast  axis  of 
the  compensator  relative  to  the  plane  of  incidence  is  X,  and  if  the  angle  between  the  fast  axis 
and  analyzer  is  7,  then  Epo/Ejo  =  1,  tan  \p  =  Rpo/Rso,  a  =  0,  and  j8  =  A 

Let  the  equation  of  the  ellipse  traced  out  by  the  electric  vector,  R,  of  the  reflected 
beam  be 


it  "2 
_ y_ 

<1* 


(18) 


where  the  components  Rx"  and  Ry"  are  referred  to  a  coordinate  axis  frame  x"y"z" 
formed  after  rotation  by  x  about  the  z'-  axis  of  the  x'-  y'-  z'  reference  frame  of  Fig.  I. 
These  components  are  related  to  the  p  and  s  components  of  k,  Rp,  and  Ri(  respectively, 
by  the  transformation  inverse  to  10: 

cos  x  sin  x  \  /  lip  \ 

09) 

-sin  x  cos  x  /  y  Ra  J 


The  analysis  of  this  elliptically  polarized  reflected  beam  is  entirely  analogous  to  that  of 
the  elliptically  polarized  beam  given  in  Section  2.1.  Thus,  the  elliptical  axes  c  and  d  can  be 
chosen  to  be 
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c  =  sin  j8 


cos2  >  sin"  x  sin  2x  cos  8 

-  (  - - +  - — 


pc. 


R  K 

p  0  so 


(20) 


and 


sin  /9  / 
/ 


sin2  x 


sin  2x  cos  fi 

R  R 

p  O  *0 


(21) 


and  the  relation  for  the  angle  of  rotation,  x,  from  the  plane  of  incidence  is 


tan  2x 


2  R  R  cos  j 0 
po  so  • 

R2  -  R2 

so  P  0 


If  Eqs.  (20)  and  (21)  are  squared,  then  added,  we  have 


sin2  $ 


o2+  d2 
c2  d2 


R2  +  R2 

so  p  ci 

R2  R2 

5  0  [tO 


(22) 


(23) 


Transformation  of  Eq.  (19)  produces 


=  Rso+«po 

The  ellipticity,  y,  of  the  ellipse  traced  out  by  the  head  of  R  is  defined  in  the  relation 

i  .  sin  v 

tan  y  =  d,  c  =  - L — 


(24) 


(25) 


Thus,  choosing  c  =  A  cos  y  gives  d  =  A  sin  7.  Of  course,  Rpo  =  R0  sin  and  Rso  =  R0  cos 
^  in  the  present  case,  and  from  Eq.  (24),  A  =  R0.  Substituting  Eqs.  (24)  and  (25)  into  Eq. 
(23)  yields 

sin"  cos^  iA  sin"  =  cos4  y  sin"  y 


or 


1"  2y  =  sin  "2i£sin2/3 


or 


sin  2 y  -  ±sin  sin  A 


(26) 
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where  the  substitution  A  =  (3  is  made.  Also,  from  Eq.  (22), 


or 


tan  2x  = 


2(RPo/R«°)C°Sfl 

1  -  U2  / R2 

P  0  /  SO 


(2  tan  lit  \ 
1  -  tan2  ff>  / 


tan  2x  =  tan  2 \j/  cos  A 


(27) 


We  next  find  the  relations  for  A  and  \j/  reciprocal  to  Eqs.  (26)  and  (27).  Squaring  Eq.  (26) 
and  noting  sin2  6  =  1-  cos2  6  produces 


cos2  2y  = 


or 

cos  2y  =  cos  2^  (1  +  ,eos2  A  tan2  2?-/ )!i 
which,  upon  substitution  of  Eq.  (27),  is 

cos  2y  =  cos  lijj  (1  +  tan^  2y)  =  cos  2|/r/'cos  2x 

or 

cos  2i}/  =  cos  2v  cos  2x 


(28) 


Eq.  (27)  may  also  be  written  as 

(sin  2x  )/  (cos  2x  )  =  j^(sin  2i^)/(cos  2t/r)  cos  aJ 

Substituting  sin  2\p  =  ±  [(sin  2y)/(sin  A)]  from  Eq.  (26)  and  for  cos  2\p  from  Eq.  (28)  gives 


sin  2x  =  + 


[(sin  2y)/(sin  A)  ^ 

cos.  2y  J 


or 


tan  A 


+ 


/  tan  2y  \ 
\  sin  2x  / 


(29) 


Corrections  for  slight  deviations  from  the  ideal  quarter-wave  compensator  phase  shift  of 
tc/2  have  been  given  by  Plumb  (Ref.  12).  We  should  note,  however  that  the  quarter-wave 
null  method  must  be  done  at  discrete  wavelengths,  which  are  determined  by  the  availability 
of  suitable  compensators.  The  next  null  method  overcomes  this  limitation. 
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2.4.3  The  Fixed  Compensator  Null  Method 

In  the  fixed  compensator  method  of  ellipsometry,  the  fast  axis  is  set  at  a  fixed  angle  to 
the  plane  of  incidence  —  usually  ±(ir/4)  —  and  the  null  occurs  by  adjusting  the  polarizer 
and  the  analyzer.  Either  the  PCSA  (Fig.  2)  or  the  PSCA  configurations  may  be  used.  The 
rather  lengthy  derivation  of  the  ellipsometric  relations  for  V'  and  A  in  the  fixed  compensator 
method  will  now  be  shown. 

For  a  PCSA  ellipsometer  arrangement,  consider  the  plane-polarized  beam  transmitted  by 
the  polarizer  with  an  electric  vector  Ej  =  EjU;  (Fig.  3)  where  Uj  js  a  unit  vector  in  the 
plane  of  the  electric  field  and  perpendicular  to  the  propagation  direction.  The  magnitude, 
Ei,  of  the  electric  field  is  the  real  part  of 


In  Fig.  3,  the  angle  P  (measured  by  the  polarizer)  is  the  angle  of  the  plane  of  the  electric  field 
with  the  plane  of  incidence.  The  beam  is  incident  upon  the  compensator,  whose  fast  axis 
makes  an  angle,  C,  with  the  plane  of  incidence. 


17 


AEDC-TR -80-51 


The  fast  and  slow  components,  Ee  and  E0,  respectively,  of  the 


A 

complex  amplitude,  E;, 


are 


Ee  =  E;-cos(C  -P)  (31) 

and 

E0  =  -E.  *  sin  (C  -  P)  (32) 

The  corresponding  amplitude  components,  EJ  and  Eg,  of  the  elliptically  polarized  beam 
transmitted  by  the  compensator  are 

E'  =  f0  •  E;  •  cos  (C  -  P)  -  eiS  (33) 


and 


-f  ■  E-  ■  sin  (C  —  P) 


(34) 


The  phase  lead,  6,  of  the  fast  component  is  attributable  to  the  compensator;  fe  and  f0 
respectively  represent  the  fractions  of  the  incident  amplitudes  of  the  fast  and  slow 
components  that  exit  the  compensator.  The  relationships  among  (1)  the  p  and  s  components 
of  the  complex  amplitude  of  the  elliptically  polarized  beam  exiting  the  compensator  and  (2) 
the  fast  and  slow  components  of  the  electric  field  are  given  by  the  transformation  in  Eq. 
(10): 


cos  C  —  sin  C  \  /  E' 


sin  C  cos  C  /  \  E ' 


iS 


E- 


fe  cos  (C  —  P)  cos  C  e  +  f0  sin  (C  —  P)  sin  C 


kfc  cos  (C  —  P)  sin  Ce  —  fQ  sin  (C  —  P)  cos  C 


(35) 


It  is  often  assumed  that  fe  =  f0;  i.e.,  the  compensator  is  perfectly  dichroic,  and  the 
differences  in  multiple  reflection  effects  along  the  two  optical  axes  are  negligible.  Deviations 
from  ideal  behavior  have  been  observed  (Ref.  13),  and  corrections  for  real  compensators  will 
be  discussed  after  the  fe  =  f0  case  is  considered. 


Let  us  assume  that  fe  =  f0  and  choose  C  =  ±(ir/4)  so  that  sin  C  =  cos  C.  From  Eq. 
(35), 


1  +  cot  Or/ 4  —  P)  c1  ® 
— 1  +  cot  Or/ 4  ~  P)  e1^ 


(36) 
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Following  Archer  and  Shank  (Ref.  14),  we  will  show  that,  for  the  ratio  tan  L  of  the 
magnitude  of  the  p  and  s  amplitudes,  i.e., 


tan  L 


-IV' 


and  for  the  phase  difference  A'  between  the  p  and  s  components,  Eq.  (36)  yields 


tan  A '  =  tan  (2P  -  tr/2)  ■  sin  8 


and 


cos  2L  =  —cos  2P  cos  8 


(37) 


(38) 

(39) 


If  the  numerator  and  denominator  of  Eq.  (36)  are  both  multiplied  by  the  complex  conjugate 
of  the  denominator  and  if  the  appropriate  trigonometric  identities  are  inserted,  we  obtain 

A  1 

col2  6  —  1  —  2i  cot  0  sin  S 
^■8  esc2  6  —  2  col  9  cos  8 


[loot2  8  -  l)2  1-  4  col 2  6  sin2  aP 
esc2  8  —  2  cot  8  cos  8 

tan  L  exp  (iA  ') 


exp 


i  tan  "1 


—2  cot  8 
cot2  9  —  1 


sin  8^ 


(40) 


where 


6  =  U/2J-P 


(41) 


From  the  identities 

tan  2x  =  (cot  2x) =  (2  cot  x)/(cot  2  x  -  1) 

and 

tan  A '  =  —tan  20  sin  8  =  tan  (—  26)  sin  8  (42) 

the  above  equation  becomes  Eq.  (38)  upon  substitution  of  Eq.  (41). 

The  moduli  of  Eq.  (40)  are  also  equal;  therefore, 

cot2  2  0+  sin2  8 


tan2  L 


(esc  2  9  —  cos  5  )2 


(43) 
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if  a  factor  of  4  cot2  0  is  extracted  from  both  the  numerator  and  denominator  of  the  squared 
modulus  and  if 


>  sin2x  =  2  siux  cos  x 

is  invoked  along  with  the  identity  in  Eq.  (42).  Noting  that 

9  .  1  -  cos  2  L 

tan'*  L  -  - 

1  +  cos  2L 

and  solving  for  cos  2L  in  terms  of  0  and  6,  we  obtain  upon  simplification 

cog  2[  _  esc2  20  -  2csc  20  cos  S  -  cos2  8  -  cot2  20  -  sin2  $ 
esc2  20  —  2csc  20  cos  5  +  cos2  S  —  cot2  20  —  sin2  8 


When  the  identities 


n  n  ?  ? 

1  =  sin.'*  x  +  cos  x  =  esc  x  —  cot  x 


are  inserted  into  Eq.  (46),  we  find 

_ _  91  cos  5  (cos  5  —  CSC  20)  . 

cos  2L  =  - — - - - —  =  -sin  20  cos  5 

esc  0  {esc  20  —  cos  5) 

Equation  (48)  becomes  identical  to  Eq .  (39)  when 

sin  20  =  cos  2P 


is  used. 


(44) 


(45) 


(46) 


(47) 

(48) 


(49) 


Now  the  relations  must  be  found  for  the  parameters  ^  and  A  in  the  fixed  compensator 
mode  of  null  ellipsometry.  Because  the  reflected  beam  is  plane-polarized,  and,  thus,  the 
phase  differences  /3,  between  the  p  and  s  components  of  the  electric  field  of  the  reflected 
beam  are  integer  multiples  of  v,  the  change  in  the  phase  difference,  A,  is  either 

A,  =  -A',  for/9  =  2Itt  (50) 


or 


A2  =  -A '+  r,  for  j8  =  (21!  -  1)  ir  (51) 

where  f  is  an  integer.  These  two  values  of  A  correspond  to  0  =  (t/4  -  P)  ranges  of 
0  <  0  s  -ir  and  -ir  <  0  <  0,  respectively,  so  there  are  two  extinction  settings  of  the 
polarizer  angle,  e.g..  Pi  and  P2.  Because  two  distinct  angles  having  the  same  tangent 
function  value  must  differ  by  an  odd  multiple  of  ir,  we  have  from  Eq.  (39), 
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P2-Pj  =  rr/2  (52) 

A  relation  for  \p  is  obtained  from  the  measurement  of  A,  the  angle  of  the  analyzer  from 
the  plane  of  incidence  at  extinction.  It  is  seen  in  Fig.  4  that  the  plane-polarized  electric 
vector,  R,  of  the  reflected  beam  is  at  the  angle  A  -  (*72),  giving  (from  Fig.  4a) 

I  R  I 


ft 


=  -cot  (Aj  -rr/2)  -  tan(-Aj) 


(53) 


for  *72  £  Ai  <  v,  and  (from  Fig.  4b) 


cot  (Ag  —  n/ 2)  =  tan  A2 


from  the  definition  of  Eq.  (15),  and  Eqs.  (37),  (53),  and  (54), 

tan  i (f  =  tan  (-A^)  cot  Lj  =  tan  A2  cot  L2 
for  the  two  sets  of  extinction  settings. 


(54) 


(55) 


\j/  can  be  related  to  measurable  parameters  from  a  consideration  of  the  orientation  of  the 
electric  vector.  In  Fig.  5,  the  electric  vector  incident  upon  the  compensator  is  pictured  for  the 
two  cases  P  =  Pi  and  P  =  P2  =  Pj  +  (rr/2).  We  can  see  that  the  fast  component  of  the 
electric  vector  for  the  case  P  =  Pj  has  an  amplitude  identical  to  the  amplitude  of  the  slow 
component  in  the  case  P  =  P2,  and  vice  versa.  Thus,  Eqs.  (53)  and  (54)  indicate 

tan  L2  =  cot  Lj  (56) 

and,  from,  Eq.  (55), 

tan ^  i/f  =  tan  (  —  A^)  tan  A2  (57) 

It  is  important  to  note  that  in  Eqs.  (33)  through  (48)  the  compensator  phase  shift,  5,  need 
not  be  rr/2.  This  means  that  during  an  ellipsometric  run,  the  wavenumber  of  the  radiation 
may  be  varied  without  replacing  the  fixed  phase  shift  compensator.  This  represents  a  distinct 
advantage  over  the  quarter-wavelength  compensator  null  method. 


If  C  were  chosen  to  be  —  (*74),  the  steps  outlined  above  would  produce  relationships  • 
among  the  ellipsometric  parameters  and  A  and  the  instrumental  settings  P  and  A,  which 
are  similar  to  Eqs.  (50),  (51),  and  (57). 
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b.  Analyzer  angle  <  tt/2 

Figure  4.  Electric  field  of  the  plane-polarized  beam 
incident  upon  the  analyzer. 
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Figure  5.  Relationship  of  the  two  possible  polarizer  angle 
values  Pi  and  P2. 

Archer  and  Shank  (Ref.  14)  have  demonstrated  empirically  that  many  compensator 
plates  have  unequal  beam  attentuations  along  the  fast  and  glow  orientations;  i.e.,  fe  *  f0.  It 
may  be  that  the  multiple  reflection  effects  within  the  plate  are  noticeably  different  along  the 
two  optical  axes  because  of  the  differences  in  the  refractive  index,  and/or  the  plates  are 
pleochroic.  This  causes  P2  -  Pi  *  ir/2,  and  it  invalidates  Eq.  (38)  for  A' .  However,  Eq.  (57) 
for  tp  is  still  correct.  Archer  and  Shank  (Ref.  14)  give  the  corrections  necessary  to  compute  A 
from  the  two  measured  null  positions  of  the  polarizer  angle,  P. 

2.4.4  Principal  Angle  Null  Methods 

The  principal  angle  of  incidence,  8,  is  defined  as  the  angle  of  incidence,  tf>ot  at  which  the 
change  of  phase  difference,  A,  between  the  p  and  s  components  is,  upon  reflection,  equal  to 
±(t/2).  It  has  been  found  that  the  nulling  schemes  discussed  above  are  usually  most 
sensitive  when  the  angle  of  incidence  is  near  the  principal  angle.  For  a  great  many  sample 
systems,  6  is  about  70  deg. 

In  most  of  the  principal  angle  of  incidence  methods,  the  angle  of  incidence  is  varied  until 
the  principal  angle  of  incidence  is  found  for  different  fixed  settings  of  the  polarizer, 
compensator,  and  analyzer  angles  (Refs.  4,  15,  16,  and  17),  Kinosita  and  Yamamoto 
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describe  one  principal  angle  of  incidence  scheme  in  which  the  compensator  is  eliminated 
(Ref.  17).  The  primary  advantage  of  these  methods  lies  in  the  ease  with  which  the  optical 
parameters  could  be  determined  for  A  =  ±  (tt/2).  However,  these  methods  require  more 
precise  measurements  of  the  angle  of  incidence  to  achieve  the  accuracy  of  other  ellipsometric 
methods. 

2.4.5  Nonextinction  Methods 

Most  milling  schemes  in  ellipsometry  require  adjustments  of  at  least  two  highly 
interdependent  instrumental  components.  Precision  measurements  often  take  considerable 
time,  even  with  the  use  of  automatic  ellipsometers.  Rapid  changes  in  surface  conditions, 
such  as  those  occurring  during  film  formation,  cannot  be  monitored.  Special  problems  may 
sometimes  be  posed  by  the  difficulty  of  obtaining  precision  compensators  and  polarizers  for 
IR  ellipsometry.  Increased  measuring  speed  and  more  accurate  IR  measurements  are  often 
obtained  through  observing  intensity  changes  in  the  ellipsometer. 

In  most  non-nulling  ellipsometric  techniques,  all  instrumental  components  except  one 
are  held  at  predetermined  settings.  The  movable  component,  usually  the  analyzer,  is  rotated 
at  a  constant  rate,  and  the  output  of  the  detector  is  electronically  monitored  for  changes  in 
intensity  and  phase.  While  changes  in  the  ellipsometric  parameters  ^  and  A  can  be  closely 
followed,  the  accuracy  of  measurement  of  the  optical  parameters  for  this  method  is  usually 
less  than  the  accuracy  attainable  by  nulling  techniques.  This  is  because  of  the  imprecision 
associated  with  measurements  of  the  varying  intensity  and  phase.  (See  Refs.  4,  8,  and  18 
through  21.) 

In  many  instances  of  nonextinction  ellipsometry,  both  the  compensator  and  its 
drawbacks  can  be  eliminated.  Such  ellipsometric  methods  are  particularly  attractive  for  IR 
work  (Refs.  8,  20,  and  21)  in  which  accurate  intensity  measurements  are  difficult  to  obtain. 
Stobie  el  at  (Ref.  20)  discuss  an  IR  ellipsometer  that,  in  some  cases,  does  not  require  such 
critical,  absolute  intensity  measurements.  Their  ellipsometer  has  a  stationary  polarizer  and 
analyzer  and  a  rotating  polarizer  (Pr0i)  1°  either  P  PFOt  S  A  or  a  P  S  Prot  A  configuration.  By 
relating  the  component  frequencies  of  the  intensity  variation  to  the  phase  shifts  of  the  p  and 
s  components  of  the  beam,  one  may  calculate  $  and  A. 

In  all  ellipsometric  intensity  measurements,  the  quadrant  in  which  a  calculated  angle  lies 
cannot  be  determined  because  only  the  absolute  values  of  amplitudes  are  observed.  This 
means  that  good  approximate  values  of  the  optical  constants  must  be  known  prior  to 
measurement  so  that  the  proper  signs  of  angles  can  be  inferred. 
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2.5  ERRORS  ASSOCIATED  WITH  ELLIPSOMETER  COMPONENTS 

The  errors  in  angular  measurement  and  orientation  and  the  spurious  effects  caused  by 
reflections  and  depolarizations  of  the  beam  can  be  ascribed  to  all  ellipsometer  components. 
Some  of  the  problems  associated  with  imperfect  compensators  are  discussed  in  Section 
2.4.3.  The  causes  and  proposed  corrections  for  many  instrumental  errors  have  been 
compiled  and  explored  by  Azzam  and  Bashara  (Ref.  4).  Error  analysis  depends  to  a  large 
extent  on  the  particular  ellipsometric  technique  used  and  the  ellipsometer  configuration. 

For  many  ellipsometric  measurements,  including  those  on  cryogenic  or  ultraclean 
surfaces,  an  environmental  chamber  surrounding  the  sample  system  is  required.  The 
chamber  windows  become  new  optical  components  in  the  optical  train  of  the  ellipsometer 
and  can  be  birefringent  under  the  stress  of  a  pressure  differential  of  one  atmosphere.  Also, 
nonplanar  and/or  nonparallel  window  surfaces  can  cause  deviations  of  the  beam  path. 

Experimentally,  the  following  procedure  may  be  used  for  vacuum  measurements  of  ^ 
and  A  (Ref.  3).  After  a  vacuum  measurement,  the  chamber  is  allowed  to  return  to  the 
ambient  room  pressure  and  temperature.  During  this  time  a  clean  surface  will  oxidize  and  a 
cryogenic  film  will  evaporate.  After  the  oxidation  or  evaporation  is  complete,  the  chamber 
may  be  pumped  out  again,  and  measurements  of  ^  and  A  repeated  for  the  new  surace  in 
vacuo.  After  the  chamber  interior  is  once  again  at  ambient  pressure  and  temperature  and  the 
windows  are  removed,  ^  and  A  are  remeasured  for  the  new  surface.  After  the  results  are 
compared  for  the  new  surface  at  the  two  experimental  conditions,  the  vacuum 
measurements  of  the  original  sample  system  can  be  corrected. 

2.6  RELATIONSHIP  OF  THE  ELLIPSOMETRIC  PARAMETERS  i  AND 
A  TO  THE  OPTICAL  PROPERTIES  OF  THE  SAMPLE  SYSTEM 

The  optical  properties  of  the  system  as  well  as  its  geometry  directly  affect  the  outcome  of 
ellipsometric  measurements  of  and  A  (that  is,  if  the  optical  properties  and  the  geometry  of 
a  sample  are  known),  the  values  of  and  A  can  be  computed.  The  inverse  problem  of 
extracting  the  sample  optical  properties  and  geometry  from  the  ellipsometric  measurements 
of  and  A  is  the  more  common  situation,  and  we  shall  consider  their  problem  for  the  two 
sample  cases  of  (1)  a  semi-infinite  ambient  and  substrate,  and  (2)  a  uniform  thin  film 
sandwiched  between  a  semi-infinite  ambient  and  a  semi-infinite  substrate.  In  what  follows 
we  shall  assume  that  all  materials  in  the  sample  system  are  isotropic.  (Anisotropic  media  are 
treated  in  Refs.  4  and  22.) 
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Figure  6  shows  the  planar  interface  of  two  semi-infinite  media  and  a  beam  incident  upon 
the  interface  at  an  angle  4>0.  Adapting  the  notation  of  Ref.  1  to  the  case  of  non-normal 
incidence  for  two  media  ns  and  nj,  we  may  write  the  complex  Fresnel  boundary  conditions 
for  the  reflection  coefficient,  rjj,  and  the  transmission  coefficient,  tjj,  as  (Ref.  23)  ■ 

"i cos  <f>j  ~  cos  <f> j 

n.  cos  A.  +  n.  cos  A. 
i  ■  J  l 


2  B .  cos 


n,  cos  <£.  +  n.  cos 

i  j  j 


rijs 


n.  cos  d>,  —  n.  cos  <A. 

t  I 

n.  cos  d>.  +  n.  cos  <£. 
1  J 


(59) 

(60) 


and 


tijs 


2  n .  cos  A. 

i  "i 


n.  cos  +  iv  cos  <f>^ 


(61) 


It  is  clear  from  Eqs.  (58)  through  (61)  that  the  p  and  s  components  do  not  behave  in  the  same 
manner  except  at  normal  incidence. 


Figure  6.  Reflection  and  refraction  for  an7 ambient-substrate  system. 
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2.6.1  The  Ambient-Substrate  System 


From  the  Fresnel  relations  [Eqs.  (58)  through  (61)],  Ditchburn  (Ref.  24)  has  derived 
equations  for  the  optical  indices  n0  and  ni  for  the  case  of  a  semi-infinite  nonabsorbing 
ambient,  i.e.,  n0  =  nQ  +  ik<,  =  nD.  These  are  (Refs.  25  and  26) 

„2  2  ■  2  , 
n  i  =  "0  s-"1  <P0 


1  -U 


tan2  cos  2ijj 


[_  (l  +  sin  2i{/'  cos  A)2 


(62) 


and 


ki  - 


no  tan  tpo  sin  tf>a  sin  2 ip  sin  A 
1  -I-  sin  2\p  cos  A 


(63) 


ip  and  A  represent  the  ellipsometric  parameters  of  a  film-free  substrate;  thus,  the  optical 
constants  can  be  found  for  an  ultraclean  substrate. 


Light  beams  in  an  absorbing  medium  such  as  a  metal  are  not  homogeneous  (Ref.  5),  so 
the  optical  constants  change  with  the  angle  of  incidence.  Vasicek  (Ref.  26)  shows  that  the 
optical  constants  at  normal  incidence  ri  =  n  +  ik,  are  preferred,  which  he  gives  as 


"l  = 


C  +  A 


(64) 


*5--^ 
1  2 


(65) 


where 


A  =  n|  —  k  |  =  n  j  -  k  ^ 

B  =  2rij  k  ^  cos  rp]  =  2nj  kj 
C2  =  A2  +  Bz 


and 


(66) 
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2.6.2  The  Ambient-Thin  Film-Substrate  System 

In  Fig.  7  it  is  assumed  that  the  thin  film  is  uniform  and  homogeneous  and  has  parallel 
boundaries  between  the  ambient  and  substrate  media.  Also,  the  substrate  and  the  ambient 
are  considered  semi-infinite  in  extent:  i.e.,  multiple  reflection  effects  within  the  substrate 
and  the  ambient  are  assumed  not  to  exist.  For  highly  absorbing  substrates  (e.g.,  metals) 
which  are  actually  thick  films,  this  assumption  seems  justified  because  the  virtual  extinction 
of  a  beam  in  the  substrate  occurs  over  pathlengths  on  the  order  of  a  wavelength  throughout 
most  of  the  IR,  visible,  and  near-UV  spectral  regions.  For  substrates  which  are  dielectric 
thick  films,  however,  multiple  reflections  at  the  dielectric  interfaces  can  occur,  and  their  ef¬ 
fects  have  been  successfully  treated  in  conventional  reflection  (Ref.  27)  and  transmission 
{Ref.  21)  spectroscopy.  Nevertheless,  we  have  yet  to  find  in  the  literature  a  treatment  of 
multiple  reflection  effects  in  a  substrate  applied  to  ellipsometric  techniques,  and  feel  that 
this  topic  may  be  an  area  of  future  interest. 

The  fundamental  relations  for  the  interpretation  of  ellipsometric  measurements  of  an 
isotropic  and  lamellate  thin  film  on  a  substrate  were  derived  by  Paul  Drude  in  1889.  The 
complex  reflection  coefficient,  roi2,  for  the  thin  film  and  substrate  in  Fig.  7  is  (Refs.  2  and 
25) 

roi2E  =  [roi£+?i2ee*p  (2ind)j/[1  +  ?oifri2f  exP  ^2iyi^J  (67) 

E  =  p,  s 


where 


Yi 


=  2i7i/n, 


^1 


(68) 


and  </>i  can  be  found  from  Eq.  (66)  if  the  ambient  is  nonabsorbing.  From  Eqs.  (13),  (14), 
and  (67),  we  may  write 


tan  tp 


e 


-i  A 


P01p  +  ?12p  e*P^y!d> 

1  +  ?01P  ?12P  exP  C2  iy  !  d) 


1  +  ?01«?12«  e*p  (2  iyt  d) 
?oi.  +  fi2S  e*p(2iyid) 


(69) 


which  is  often  called  the  “exact”  Drude  equation.  Thus,  \p  and  A  are  complicated  functions 
of  nine  parameters:  v,  <f>OJ  and  d,  and  the  real  and  imaginary  parts  of  n0,  rij ,  and  n2. 
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Figure  7.  Reflection  and  refraction  for  an  ambient-thin  film- 
substrate  system. 


The  inversion  of  the  real  and  imaginary  parts  of  Eq.  (69)  to  solve  for  any  of  the  nine 
parameters  cannot  be  done  precisely.  Laborious  iterative  algorithms  must  be  used  for  the 
inversion,  calculations  impractical  until  the  advent  of  modern  digital  computers.  For  this 
reason,  Drude  and,  later,  Tronstad  developed  approximate  expressions  for  the  quantities 
A-A  and  vW1  which  were  linear  in  film  thickness,  d,  but  accurate  only  for  very  thin  films 
(d  <  <  X).  Here,  the  bar  indicates  the  parameters  for  the  film-free  case  of  Section  2.6. 1 .  For 
certain  instances,  Burge  and  Bennett  (Ref.  29)  report  large  discrepancies  between  the  optical 
parameter  values  determined  from  the  approximate  relations  and  those  determined  from  Eq. 
(69).  Most  present-day  investigators  make  use  of  computer  algorithms  for  obtaining  the 
desired  parameters  from  Eq.  (69).  The  algorithms  vary  according  to  the  parameters  sought, 
the  information  known,  and  the  criteria  for  convergence. 

In  principle,  all  nine  parameters  can  be  determined  if  at  least  nine  independent 
ellipsometric  measurements  are  made.  However,  in  nearly  all  cases,  the  parameters  v,  <£o, 
and  nQ,  the  ambient  refractive  index,  are  known.  Of  the  remaining  five  parameters,  the 
substrate  complex  refractive  index  n2  =  n2  +  ik2  is  often  also  obtained  from  previous 
measurements,  leaving  just  the  three  film  parameters  nlf  kj,  and  thickness,  d,  to  be 
determined.  In  a  typical  ellipsometric  experiment,  only  two  quantities  are  measured,  a  single 
value  of  ^  and  one  of  A.  Thus,  unless  the  film  is  known  to  be  nonabsorbing  or  unless  its 
thickness  is  known,  the  three  unknown  film  parameters  cannot  be  found  from  a  single 
typical  ellipsometer  run,  and  additional  information  must  be  sought. 
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2. 6. 2.1  Case  of  Only  Two  Unknown  Parameters 

If  only  two  of  the  original  nine  parameters  in  the  ambient-thin  film-substrate  system  are 
unknown,  the  unknown  parameters  may  be  determined  from  a  measured  value  and  a 
measured  A  value.  In  principle,  any  two  of  the  nine  parameters  may  be  unknown,  but 
frequently  the  two  unknowns  are:  (l)  ni  and  d  of  a  transparent  film,  (2)  ni  and  k|  of  an 
absorbing  film,  or  (3)  the  optical  constants  r\2  and  ki  of  the  substrate. 

Several  authors  (Refs.  30  and  31)  have  written  algorithms  appropriate  for  cases  1  and  2 
above.  McCrackin  (Ref.  30)  describes  a  FORTRAN  program  that  will  compute  the  two 
unknown  parameters,  either  nj  or  k|,  or  ni  and  d|,  from  the  ellipsometer  settings  of  the 
polarizer  and  analyzer  angles,  P  and  A.  The  program  corrects  for  differences  in  the 
attenuation  of  the  fast  and  slow  components  of  a  beam  passing  through  an  imperfect 
compensator  (Section  2.4.3),  and  misalignment  of  the  sample  surfaces. 

2.6.2.2  Case  of  Three  or  More  Unknown  Parameters 

If  three  or  more  parameters  (e.g.,  kj,  and  d  of  an  absorbing  film)  are  to  be 
determined  from  ellipsometric  measurements,  it  may  be  possible  to  obtain  sufficient 
independent  observations  of  ^  and  A  by  varying  either  (1)  the  film  thickness,  d,  (2)  the 
ambient  medium,  nG,  (3)  the  substrate  medium,  rh,  (4)  the  angle  of  incidence,  <fr0  (multiple 
angle  of  incidence  (MAI)  eliipsometry),  (5)  the  film  material,  nj,  or  (6)  any  given 
combination  of  d,  n0,  n i ,  n2,  and  McCrackin  and  Colson  (Ref.  32)  have  studied 
applications  of  schemes  1  through  4  to  find  the  three  absorbing  film  parameters  m,  kj,  and 
d.  They  concluded  that  because  of  experimental  uncertainties  the  variation  of  film  thickness 
method  —  method  1  —  works  only  for  film  thicknesses  greater  than  about  0.010  pm,  even  if 
ni  and  k|  are  assumed  to  not  change  with  thickness.  They  found  that  the  methods  employing 
different  ambients  or  substrates  (methods  2  and  3)  were  most  likely  to  yield  unique  film 
parameters.  Of  these  two  methods,  method  2  —  the  variation  of  the  ambient  (usually  a 
fluid)  —  would  ordinarily  be  easier  to  apply  than  method  3  because  of  the  difficulty 
inherent  to  method  3  of  depositing  identical  films  on  different  substrates. 

Several  authors  have  discussed  multiple  angle  of  incidence  (MAI)  eliipsometry,  method  4 
above  (Ref.  33).  Azzam  and  Bashara  (Ref.  4)  provide  general  criteria  for  choosing  the  angles 
of  incidence  and  the  least-squares  converging  function  to  be  minimized  if  the  number  of 
independently  measured  quantities  exceeds  the  number  of  unknown  parameters.  For  a 
simulated  SiO  film  on  a  Si  substrate  in  which  five  parameters  (nj ,  kj,  d,  n2,  and  k2)  were 
found  by  the  least-squares  algorithm  of  Marquardt  (Ref.  34),  Azzam  and  Bashara  noted 
more  rapid  convergence  for  very  thin  films  ( <  0.010  jim)  than  for  thicker  films  ( >  0. 10  ftm). 
They  also  concluded  that  for  all  film  thicknesses,  convergence  to  the  proper  parameter 
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values  failed  unless  fairly  accurate  (±15  percent)  initial  estimates  of  the  parameter  values 
were  used. 

i 

To  date,  we  have  not  seen  any  analyses  of  ellipsometric  schemes  in  which  a  given 
combination  of  the  parameters  4>at  d,  n0,  ni,  and  n2  is  allowed  to  vary.  In  principle,  at  least, 
such  schemes  utilizing  Marquardt’s  algorithm  (Ref.  34)  are  possible  and  show  promise  for 
future  applications. 

Additional  information  may  come  from  non-ellipsometric  methods.  For  example,  in  the 
case  of  1R  ellipsometric  measurements  of  an  unknown  absorbing  thin  film,  the  method  of 
Templemeyer  and  Mills  (Ref.  35)  may  be  used  for  an  independent  measurement  of  film 
thickness,  and  the  parameters  n]  and  ki  may  be  found  from  a  single  measurement  of  $  and 
one  of  A. 

2.7  COMMERCIALLY  AVAILABLE  ELLIPSOMETRIC  INSTRUMENTATION 

The  two  major  suppliers  of  ellipsometric  equipment  are  Rudolph  Research  and  Gaertner 
Scientific  Corporation.  Both  manufacturers  offer  the  two  styles  of  ellipsometers  listed  in 
Table  1:  (1)  production  models  that  can  be  fully  automated  and  (2)  research  instruments  that 
can  be  used  to  make  more  precise  measurements  and  that  can  be  partially  automated  with 
optional  components.  The  accessories  and  warranties  offered  by  both  firms  are  essentially 
the  same.  All  commercially  available  ellipsometers  are  in  the  PCSA  configuration. 

Table  1.  Commercially  Available  Ellipsometers 


Ellipsometer 

Type 

Gaertner 

Rudolph  Research 

Production 

Manual 

L117 

RR.100 

Automatic 

L116 

AutoEL-II  and 
AutoEL-III 

Research 

L118GT  (0.1-deg 

436  (Vertical 

Precision) 

Sample) 

L119  (0.01-deg 

437  (Horizontal 

Precision) 

Sample) 
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2.7.1  Commercially  Produced  Ellipsometers 

The  production  ellipsometers  are  designed  for  easy  and  fast  measurements  of  the 
thickness,  d,  and  refractive  index,  m,  of  a  thin  transparent  film  on  a  known  substrate,  and 
at  the  He-Ne  laser  wavelength  of  0.6328  jim.  The  two  or  three  allowed  angles  of  incidence 
are  usually  pin  adjustable.  The  automatic  models  achieve  rapid  response  either  with  a 
continuously  rotating  analyzer  or  with  stepper-motor-driven  polarizers  and  analyzers 
(Ref.  36).  This  allows  the  use  of  a  programmable  calculator  to  compute  $  and  A,  and,  often, 
ni  30,1  d.  The  production  ellipsometers  are  not  designed  For  measurements  of  sample 
systems  enclosed  in  environmental  chambers  but  can  accommodate  liquid  sample  systems  on 
their  standard  horizontal  mounts. 

Both  manufacturers  supply  optional  optics  that  can  reduce  the  beam  cross-sectional 
diameter  to  26  /an.  This  allows  a  high-resolution  mapping  of  film  thickness  across  a 
horizontally  mounted  substrate.  Indeed,  the  Gaertner  LI  15  system  includes  an  LI  16 
ellipsometer  and  associated  electronics  that  can  automatically  provide  a  hard  copy  mapping 
of  thin-film  thickness  over  a  predetermined  area  of  a  substrate. 

The  Gaertner  LI  16  automatic  production  ellipsometer  is  essentially  an  LI  17  manual 
production  ellipsometer  with  a  rotating  analyzer.  (See  Section  2.4.5.)  The  angles  of 
incidence  can  be  30,  50,  or  70  degrees.  A  Hewlett-Packard  HP9825  calculator  can  process 
the  ellipsometric  measurements  to  find  the  index  of  refraction,  m,  and  the  thickness,  d,  of  a 
nonabsorbing  film  on  a  known  substrate.  The  calculator  can  also  control  optional  stepper- 
motor-driven  scanning  stages  for  the  sample  that  allow  different  sampling  points  anywhere 
within  a  5-in.-diam  area  on  the  sample.  Provisions  have  just  become  available  for  adding  the 
new  programs  to  the  calculator  for  multiple  transparent  films  and  for  absorbing  films.  In 
addition  to  having  the  LI  16  with  HP9825  and  the  motor-driven  scanning  stages,  the  LI  15 
system  also  has  a  Hewlett-Packard  HP7225A  plotter  to  permit  a  graphical  display  of 
thickness  scans.  The  optics  for  beam  size  reduction  must  be  factory  installed  on  all  Gaertner 
production  models  and  are  not  removable. 

The  AutoEL  series  of  ellipsometers  given  in  Table  1  comprises  the  automatic  production 
models  built  by  Rudolph  Research.  A  manual  production  model,  RR100,  is  also  built  by 
Rudolph  Research,  but  it  is  less  accurate.  The  AutoEL  ellipsometers  are  null-seeking 
(Section  2.4.3)  and  use  a  microcomputer  to  drive  stepper  motors  to  select  the  proper 
polarizer  and  analyzer  settings.  The  AutoEL-lII  ellipsometer  is  a  more  advanced  version  of 
the  lower-priced  AutoEL-II,  but  both  achieve  the  measuring  accuracies  of  the  Gaertner 
LI  16  and  LI  17  models.  The  AutoEL-II  is  designed  primarily  for  automatic  measurement  of 
the  index  of  refraction,  ni,  and  the  thickness,  d,  of  a  thin  transparent  film  on  a  known 
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substrate.  The  AutoEL-IlI  has  the  added  capabilities  of  an  accessible  computer  memory  and 
can  find  the  thicknesses  and  the  refractive  indices  of  a  double  layer  of  nonabsorbing  films  or 
the  optical  constants  of  bare  absorbing  substrates.  Other  computer  programs  may  be  added 
to  the  AutoEL  series. 

In  the  Rudolph  Research  AutoEL  series  of  ellipsometers,  the  optional  translatable 
sample  stages  are  manually  adjustable,  as  opposed  to  the  (optional)  motor-driven  stages  for 
the  Gaertner  LI  16  ellipsometer.  However,  the  beam  reduction  optics  of  the  AutoEL  series 
are  removable,  and  the  angles  of  incidence  may  be  pin-locatable  or  set  with  alignment 
prisms.  The  customer  may  specify  any  set  of  two  or  three  incident  angles  between  60  and  90 
deg  for  either  pin-set  or  prism-set  angles. 

2.7.2  Commercial  Research  Ellipsometers 

The  research  ellipsometers  from  Gaertner  are  the  L118GT  and  LI  19  models  (Table  1). 
These  are  basically  modified  spectrometers  (Ref.  13)  with  vertical  substrate  mounts  which 
can  be  translated  as  well  as  rotated.  An  optional  support  for  the  entire  ellipsometer  to 
accommodate  horizontal  samples  (e.g.,  liquid)  is  available;  however,  the  ellipsometer  must 
be  modified  at  the  factory.  Component  settings  can  be  read  to  within  0. 1  deg  on  the  LI  18GT 
and  to  within  0.01  deg  on  the  LI  19.  The  minimum  angle  of  incidence  of  both  models  is  22.5 
deg.  The  standard  ellipsometer  components  are  manually  adjusted,  but  an  optional  rotating 
analyzer  is  available  for  either  model.  Light  sources  (e.g.,  He-Ne  lasers,  or  mercury  lamps 
with  filters),  a  Babinet-Soleil  compensator,  and  detectors  are  also  optional  items. 

The  LI  19X  ellipsometer  is  the  large  based  version  of  the  LI  19  which  will  allow  the  use  of 
environmental  chambers  surrounding  the  sample,  which  must  be  mounted  vertically. 
Factory-installed  beam  reduction  optics  (to  26  yxm  diameter)  are  available  only  for  the 
LI  19X.  Extended  spectral  range  measurements  (from  0.215  to  3  ^m)  can  be  made  with  the 
L1I9XUV,  an  L119X  model  with  non-glass  optics. 

The  436  research  ellipsometer  for  vertically  mounted  samples  from  Rudolph  Research 
has  the  measuring  precision  of  the  Gaertner  LI  19.  The  Rudolph  Research  437  model  is  a 
sturdier  version  of  the  436  ellipsometer.  The  437  is  designed  for  horizontally  mounted 
samples.  The  models  are  identical  in  all  respects  other  than  ranges  of  the  angle  of  incidence, 
which  are  25  to  90  degrees  for  the  436  and  50  to  90  degrees  for  the  437.  Because  nearly  all 
ellipsometric  measurements  are  made  at  angles  of  incidence  near  the  principal  angle,  G, 
which  is  usually  about  70  deg  (Section  2.4.4),  the  more  restricted  incidence  angle  range  of 
the  437  would  present  a  problen  only  rarely,  even  in  multiple  angle  of  incidence  (MAI) 
ellipsometry  (Ref.  4).  A  mercury  light  source,  filters,  apertures,  sample  holder,  and  a 
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photomultiplier  tube  (PMT)  light  detection  system  are  supplied  as  standard  equipment  on 
the  436  and  437.  These  items  are  optional  for  the  Gaertner  research  ellipsometers. 

The  optional  items  for  the  Rudolph  Research  research  ellipsometers  include  He-Ne  laser 
and  stabilized  tungsten-iodine  sources,  Babinet-Soleil  compensators,  rotating  analyzers,  and 
digital  readout  rotators  for  the  polarizer  and  analyzer.  The  rotating  analyzer  provides  faster 
measurements  but  is  less  accurate  and  requires  more  frequent  maintenance  than  the 
manually  adjusted  analyzers.  Environmental  chambers  can  be  used  with  the  436  and  437  if 
the  optical  components  of  the  ellipsometer  are  mounted  on  optional  optical  benches;  thus  it 
is  not  necessary  to  buy  a  special  model  for  environmental  chambers,  which  is  the  case  for  the 
Gaertner  line  of  research  ellipsometers.  Beam  size  reduction  optics  are  also  available  that 
can  produce  a  beam  diameter  of  2.5  jtm,  one-tenth  the  size  of  the  smallest  diameter  available 
for  the  corresponding  Gaertner  optics. 

2.7.3  Recommendations  Concerning  Ellipsometric  Equipment 

The  limitations  of  commercial  production  ellipsometers  must  be  considered  when 
choosing  between  a  production  ellipsometer  and  a  research  ellipsometer.  These  limitations 
are 

1 .  The  production  instruments  are  less  accurate  but  faster  and  easier  to  use. 

2.  All  measurements  must  be  at  0.6328  /j.m  and  are  restricted  to  a  few  fixed  angles  of 
incidence. 

3.  There  are  more  restrictions  in  the  choice  of  sample  systems  that  can  be  analyzed 
by  production  instruments,  but  complete  data  analysis  systems  are  available  for 
them. 

4.  No  environmental  chamber  measurements  can  be  made  with  production 
ellipsometers. 

If  the  AEDC  chooses  a  production  ellipsometer,  we  recommend  that  the  Gaertner  LI  17 
manual  production  ellipsometer  be  selected.  We  base  our  choice  partly  on  the  assumption 
that  speed  of  measurement  is  not  important.  The  LI  17  can  be  factory-fitted  with  beam 
reduction  optics  and  sample-positioning  devices;  also,  all  of  the  data  reduction  capabilities 
of  the  Gaertner  LI  16  automatic  production  ellipsometer  are  available. 

However,  if  an  automatic  production  ellipsometer  is  desired  at  the  outset,  we 
recommend  the  Rudolph  Research  AutoEL-III  over  the  comparably  priced  Gaertner  LI  16 
system.  The  flexibility  of  the  AutoEL-III  is  slightly  greater  than  that  of  the  Gaertner  LI  16 
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(e.g.,  the  AutoEL-IlI  has  user-removable  beam  reduction  optics)  which  would  probably 
offset  the  faster  data  analysis  time  claimed  by  Gaertner  for  the  LI  16.  [The  analysis  times  for 
a  single  film  are  given  as  four  seconds  for  the  LI  16  (Ref.  37)  and  15  seconds  for  the  AutoEL- 
III  (Ref.  38).]  We  would  recommend  the  much  more  expensive  Gaertner  LI  15  system  only 
for  the  case  requiring  profiles  of  film  thickness  on  a  substrate. 

The  authors  strongly  recommend  the  Rudolph  Research  436  or  437  research  instruments 
if  the  AEDC  is  to  acquire  a  research  ellipsometer.  In  fact,  unless  the  AEDC  has  firm  plans  to 
never  measure  liquid  samples,  we  feel  that  a  437  should  be  purchased  because  of  its  sturdier 
construction. 

While  the  Gaertner  LI  19  research  ellipsometer  offers  the  same  capabilities  as  either  the 
Rudolph  Research  436  or  437,  changes  in  ellipsometer  configurations  are  more  awkward  and 
time-consuming  for  the  LI  19  model.  For  instance,  versions  of  the  LI  19  that  can  measure 
horizontally  mounted  samples  or  those  that  can  accommodate  environmental  sample 
chambers  must  be  assembled  at  the  factory.  In  contrast,  the  437  is  designed  especially  for 
horizontal  samples  (of  both  liquids  and  solids),  and  either  the  436  or  the  437  can  be 
converted  to  vacuum  instruments  by  the  user,  with  use  of  optional  optical  benches. 
Furthermore,  there  is  no  stock  version  of  the  LI  19  available  that  can  rival  the  versatility  of 
the  437  with  its  optical  benches  that  can  measure  horizontally  mounted  samples  in  an 
environmental  chamber. 

2.8  SUMMARY  ON  ELLIPSOMETRY 

In  the  near-IR  to  near-UV  spectral  region,  ellipsometric  methods  are  comparable,  or 
superior,  to  conventional  transmission  and  reflection  spectroscopic  techniques  for  obtaining 
both  the  optical  constants  of  thin  films  and  substrates  and  the  thicknesses  of  films.  The  wide 
variety  of  ellipsometric  instrument  configurations  creates  a  myriad  of  analysis  schemes,  not 
all  of  which  have  been  explored.  The  possibilities  of  making  ellipsometric  measurements  in 
the  IR  are  being  explored  by  a  number  of  researchers;  these  measurements  could  be  made  if 
existing  commercial  ellipsometers  were  modified. 

Only  PSCA  ellipsometers  are  commercially  available,  and  these  are  of  two  types, 
production  and  research.  The  less  accurate  production  ellipsometers  can  be  fully  automated 
to  determine  the  optical  parameters  of  simple  systems,  but  only  at  a  wavelength  of  0.6328  pun. 
With  optional  optics,  the  manual  research  ellipsometers  can  operate  in  the  spectral  region 
from  0.210  to  about  3  /xm.  The  commercial  ellipsometers  rated  best  by  the  authors  are:  (1) 
the  Gaertner  LI  17  manual  production  ellipsometer,  (2)  the  Rudolph  Research  AutoEL-lII 
automatic  production  ellipsometer,  and  (3)  the  Rudolph  Research  437  research  ellipsometer. 
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3.0  INTERNAL  REFLECTION  SPECTROSCOPY  (IRS)  OF  FILMS 

In  most  ellipsometric  measurements  (Section  2.0)  or  conventional  reflectance 
experiments,  a  refracted  ray  moving  away  from  an  interface  of  two  media  removes  part  of 
the  power  of  the  incident  beam  that  would  otherwise  go  to  the  reflected  beam.  However,  if 
the  two  media  are  transparent  and  the  incident  beam  is  in  the  medium  of  the  higher 
refractive  index,  then  a  beam  incident  at  angles  greater  than  a  specific  “critical”  angle  will 
not  cause  a  refracted  beam  to  appear  in  the  lower  index  medium,  and  the  incident  and 
reflected  beams  will  have  equal  power.  This  well-known  phenomenon  is  termed  total 
internal  reflection  (T1R)  or  simply  total  reflection.  The  word  “internal”  is  usually  applied 
because  total  reflection  is  often  seen  when  the  incident  beam  is  in  glass  or  water  and  the 
incident  medium  is  surrounded  by  or  has  a  boundary  with  the  ambient  air. 

Issac  Newton  noticed  that  when  a  beam  undergoing  total  reflection  is  within  a  glass 
prism,  a  second  prism  brought  sufficiently  close  to  the  first  (but  not  in  contact)  will  cause 
propagation  of  light  within  the  second  prism  and  a  corresponding  reduction  in  the  reflected 
beam  intensity  in  the  first  prism.  This  is  known  as  frustrated  total  reflection  (FTR)  energy 
coupling,  and  is  useful  in  the  precise  measurement  of  small  distances  (Ref.  39).  FTR 
coupling  clearly  indicates  that  an  electromagnetic  field  (the  “evanescent”  field  or  wave) 
exists  in  the  rarer  (lower  refractive  index)  medium  even  under  the  conditions  of  total 
reflection.  However,  because  the  intensities  of  the  totally  reflected  beam  and  the  incident 
beam  are  equal,  it  must  be  true  that  the  time  average  energy  of  the  evanescent  wave  must  be 
zero.  (See  Section  3.1.) 

Another  method  for  “frustrating”  total  reflection  is  to  replace  the  rarer  transparent 
medium  with  a  rarer  absorbing  medium.  Here,  the  dipole  interaction  of  the  evanescent  wave 
with  the  molecules  of  the  rarer  medium  causes  a  loss  in  the  power  of  the  reflected  beam.  This 
energy  coupling  process,  called  attenuated  total  reflection  (ATR),  reveals  molecular 
absorption  spectra  very  similar  to  transmission  spectra,  and  forms  the  basis  for  internal 
reflection  spectroscopy  (IRS). 

The  research  of  Fahrenfort  (Ref.  40)  and  Harrick  (Ref.  41)  began  the  present  interest  in 
IRS  methods.  Their  work  showed  that  IRS  was  more  sensitive  than  conventional  reflection 
methods,  especially  for  weaker  absorptions.  Harrick’s  text  (Ref.  42)  describes  the  theory  of 
IRS  and  its  numerous  applications,  and  has  become  the  standard  authority  on  internal 
reflection  spectroscopy. 

Unlike  ellipsometry,  IRS  techniques  can  be  readily  applied  to  measurements  throughout 
the  IR  spectral  region.  Other  attributes  include  (Ref.  43):  (1)  enhanced  sensitivity  to  the 
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optical  constants  of  materials  over  other  reflectance  methods;  (2)  the  ability  to  selectively 
test  the  optical  thickness  of  very  thin  films  and  surface  states  because  of  the  small 
penetration  depth  of  the  evanescent  wave;  (3)  the  low  light  loss  per  reflection  for 
nonabsorbing  samples,  thereby  making  multiple  reflection  techniques  feasible  even  for  very 
thin  films;  and  (4)  the  determination  of  the  orientations  of  molecules  at  a  surface. 

V 

3.1  THEORY  OF  TOTAL  REFLECTION  AT  THE  INTERFACE 

OF  TWO  SEMI-INFINITE  MEDIA 

From  Maxwell's  laws  of  electromagnetism  one  may  show  that  at  the  interface  of  two 
media,  the  normal  components  of  both  the  magnetic  induction  field,  B,  and  the  electric 
displacement  field,  D  =  ft2  E,  are  the  same  on  either  side  of  the  boundary  (Ref.  25).  The 
Fresnel  relations,  Eqs.  (58)  through  (61),  follow  from  these  boundary  conditions  as  well  as 
from  Snell’s  law. 

In  general,  Snell’s  law  concerning  the  angles  made  by  the  incident  and  refracted  beams  to 
the  normal  of  an  interface  of  two  media  of  complex  indices  of  refraction,  n;  and  rij,  may  be 
written 


n;  sin  <f>j  =  nj  sin  (70) 

For  two  absorbing  media  the  complex  arguments  $j  and  ‘tj  have  complicated  relationships 
with  the  angles  of  incidence  and  refraction.  Even  if  the  two  media  are  transparent,  the 
argument  corresponding  to  the  smaller  (real)  index  of  refraction  must  be  complex  for 
conditions  in  which  its  sine  function  value  is  greater  than  unity. 

First  the  simpler  case  of  total  reflection  at  the  interface  of  two  semi-infinite  transparent 
media  will  be  considered.  Then  the  effects  caused  by  replacing  the  rarer  transparent  medium 
with  a  rarer  absorbing  one  will  be  discussed. 

3.1.1  Total  Reflection  for  Two  Transparent  Media 

i 

For  the  two  semi-infinite  transparent  media  shown  in  Fig.  8.  Snell’s  law  [Eq.  (70)1  is 

rij  sin  =  02  sin  $2  (71) 

A 

The  argument  4*2  is  real  (and  equals  <62  >  <i>)  only  if  <t>  1  ^  <£120  the  critical  angle  for 
transparent  media  1  and  2.  The  critical  angle,  d>  12c,  is  the  value  of  4>i  corresponding  to 
<t>2  =  ir/2,  so 

<p]2c  =  sin" J  (□  2/n  j  )  =  sin"1  n21  (72) 
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where  n2i  is  the  relative  refractive  index  for  media  1  and  2  and  is  defined  as 

n21  =  11 2 /n  1  (73) 

If  <f>\  >  <j>  J2c.  then  <£2  has  no  physical  significance,  and  total  reflection  occurs  at  the  interface 
if  the  incident  beam  is  in  the  denser  medium,  nj.  (Note  that  if  the  incident  beam  is  in  the 
rarer  medium,  n2,  then  <£2  is  always  less  than  or  equal  to  ir/2  and  <j>\  can  never  be  greater  than 
4>12c-) 

3. 1.1.1  Representation  of  the  Electric  Fields 


The  electricfields  of  the  incident,  reflected,  and  refracted  beams  are  represented  in  Fig.  8 
by  the  vectors  E,  R,  and  T,  respectively.  From  the  orientations  of  the  vectors  pictured  in  Fig. 
8  and  adopting  the  notation  appearing  in  Eqs.  (1),  (2),  (3),  (11),  and  (12),  we  may  write 


and 


E  =  Re 


f  T  /•*  -  \  -  iet  -  ia. 

I  L(U*COS*l-U;Sin,MEpo«  P+U«  e 


y  so 


] 

•  exp  ico  ^  t - -  (x  sin  0^  +  z  cos  j  )  J  j1,  J 

R  =  Re  {[_(“x  cos  4>1  +  5zsin  «0Rp°  e  ^  +  “y  Rso  e  ^  ] 

i  ia>  -  — -  (x  sin  —  z  cos  j1  ^ 


•  exp 


T  =  Re  j  [(3,  cos  S>2  -Sa  sin  4>2)  Tpo  e  ^  +  uy  Tso  e 
•  exp^-ityj^t--^  (x  sin  4>2  +  z  cos  J  ^  ^ 


(74) 


(75) 


(76) 


In  the  above  equations,  U]  (for  1  =  x,  y,  or  z)  refers  to  the  unit  vectors  along  the 
corresponding  axis  directions  pictured  in  Fig.  8.  That  is,  the  z-axis  is  directed  normally  into 
the  rarer  medium,  n2,  at  the  point  of  incidence  on  the  interface,  and  the  x-axis  lies  along  the 
interface  in  the  plane  of  incidence.  $2  is  real  (and  equals  4>i)  only  if  <t>\  is  less  than  the  critical 
angle,  </>12c. 
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Figure  8.  Parallel  and  perpendicular  components  of  the  electric  field 
of  the  incident  (E),  reflected  (R),  and  refracted  (T)  beams 
at  the  interface  of  two  transparent  media  when  the 
incident  beam  is  in  the  optically  denser  medium. 

3. 1.1. 2  The  Fresnel  Relations 


The  amplitudes  Epo,  Eso,  Rpo,  RS0)  Tpo,  and  Tso  in  Eqs.  (74),  (75),  and  (76)  are  real  and 
constant  under  steady-state  conditions.  It  is  convenient  to  define  corresponding  constant 
complex  amplitudes  as  follows: 


-  ia 

p 
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The  Fresnel  relations  are  obtained  by  matching  electric  field  components  on  either  side 
of  the  interface  at  the  point  of  incidence  (x  =  y  =  z  =  o).  From  Eqs.  (74)  through  (77),  we 
see  that  the  transmission  and  reflection  coefficients  are  simply  the  ratios  of  the  complex 
amplitudes  because  the  time  dependences  in  Eqs.  (74)  through  (76)  are  identical  at  the  origin 
of  the  x-y-z  coordinate  system,  the  point  of  incidence!  Thus,  when  Snell’s  law,  Eq.  (72),  and 
Eq.  (77)  are  put  into  Fresnel’s  relations,  [Eqs.  (58)  through  (61)1,  the  modified  Fresnel 


relationships  result: 

“II 

tan  (d>  j 

-•a) 

s  in  tji| 

cos  tj> j  -  sin  $2  ct>s  ®2 

(78) 

r‘2p '  *i ' 

tan 

+  $2) 

sin 

A  A 

cos  <f> J  +  sin  4>2  cos 

«1  * 

-sin  (<£, 

A 

-•a) 

sin 

A 

cos  4>z  -  cos  (Aj  sin  -$2 

(79) 

V 

Fl2s  =  17  " 

sin 

sin  j 

A  n 

cos  $2  -  cos  0^  sin  4>2 

2cos 

sin  4>2 

2cos^j  s in  $2 

/o  f\\ 

"  a.  ^  ^ 

Ejj  sin  {0L  +  )  cos  (0^  —  Og)  sin  cAj  cos  sin<t„  cos 


and 


2  c os  r A  ^  sin  $2 
sin  (<£j  +  <&2  ) 


2  cos  sin  4>„ 


sin  cAj  cos  $2  +  sin  4*2  cos 


(81) 


The  Fresnel  relations 
Wolf  (Ref.  25). 


involving  sums  or  differences  of 


4>]  and  4>2  are  given  by  Born  and 


3. 1,1.3  The  Energy  Conditions  of  Total  Reflection 

Total  reflection  of  a  beam  can  be  shown  to  occur  at  the  interface  of  two  transparent 
media  when  the  angle  of  incidence,  <£i,  exceeds  the  critical  angle,  <£i2cJ  i.e. ,  all  the  power  of 
the  incident  beam  appears  in  the  reflected  beam  when  >  012c. 

When  <f>i  >  <Ai2c>  Eq.  (72)  gives 

A 

sin  $2  =  8>n  0l/n21  >  *  (82) 


Thus, 


COS 


4*2  “  £  1  —  8*n2  4*2]  =  ‘  [  S'n2  -  ^  ] 

-KM' 


(83) 
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The  positive  square  root  in  Eq.  (83)  must  be  chosen  so  that  the  electric  field  will  not  appear 
to  increase  exponentially  with  z  in  the  rarer  medium.  Substituting  Eqs.  (82)  and  (83)  into  the 
last  relation  presented  in  each  of  Eqs.  (78)  through  (81),  produces,  for  <f>i  >  <f> |2C. 


-r12p  “ 


'12s 


(84) 


(85) 


and 


2  iij  j  cos  ^ 


n^cos  <f>l  +  i  |js in2 


2  cos  <j> i 


‘12s 


4> i  +  ‘  [si"2  <f>i  ~  n2l] 


(86) 


(87) 


The  magnitudes  of  the  denominators  and  numerators  of  both  Eqs.  (84)  and  (85)  are 
identical,  giving 


and 


(88) 

(89) 


That  is,  all  of  the  intensity  of  the  parallel  component  of  the  incident  beam  appears  in  the 
parallel  component  of  the  reflected  beam  when  <t>  >  ^i2C,  and  likewise  for  the  perpendicular 
components. 


3. 1.1. 4  Phase  Shifts  Upon  Reflection 


For  total  reflection,  the  expressions  for  ri2p  and  ri2s  can  be  written  [Eqs.  (77),  (78),  (79), 
(88),  and  (89)]  as 

-i(j8g-Og) 


r12f 


rI  2  0 


i  f  =  p. 


(90) 
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where 


8l  =  Pi~at 


(91) 


The  form  of  r12f  is  shown  in  Eqs.  (84),  (85),  (90),  and  (91)  to  be 


where 


ri2 £  - 


ag  -  i  b 
ag  +  i  b 


e 


-2ig£ 


Sg  bji 

tan  a  n  =  tan  — —  =  -  ;  ^  =  P,  s 

1  2  a£ 


The  combination  of  Eqs.  (84)  and  (85)  with  (92)  and  (93)  yields 


(92) 


(93) 


and 


n 


2 

21 


COS 


COS  cb-j 


(94) 


(95) 


Therefore,  the  changes  in  polarization  at  total  reflection  are  not  the  same  for  each 
component,  and,  in  general,  the  reflected  beam  is  elliptically  polarized  even  if  the  incident 
beam  is  plane  polarized. 


It  is  interesting  to  see  what  the  ellipsometric  parameters  ^  and  A  =  <5p  -  6S  become  under 
the  conditions  of  total  reflection  for  two  semi-infinite  transparent  media  (Section  2.0).  From 
Eq.  (11), 


P 


rl  =  P  =  e-i  A 

r12s 


which  upon  comparison  with  Eqs.  (12),  (94),  and  (95)  gives 


(96) 


lan  i^=l 


(97) 
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and 


tan  —  —  tan  — 

2 _ 2_ 

8  8, 

P  8 

]  tan  —  tan  — 


^[sin2  n2  , ]  ' 
sin2 


(98) 


So,  for  ellipsometric  measurements  of  two  semi-infinite  transparent  media  under  the 
conditions  of  total  reflection,  $  is  already  determined  from  Eq.  (97)  to  be  tt/4  plus  an  integer 
multiple  of  r,  and  only  A  needs  to  be  found.  The  relative  index  can  be  calculated  from. 
Eq.  (98). 


3.1. 1.5  The  Nature  of  the  Evanescent  Field 

One  may  evaluate  Eqs.  (86)  and  (87)  for  it^pl2  and  |'ti2S  [2  and  see  that  neither  quantity 
is  zero,  which  indicates  that  some  power  must  enter  the  rarer  medium,  forming  the 
evanescent  field.  If  the  time  average  power  in  the  rarer  medium  is  computed,  however,  one 
finds  it  is  zero  (Ref.  25),  and  the  evanescent  field  has  the  form  of  a  standing  wave  (Ref.  46). 

A  A 

The  form  of  the  evanescent  wave  is  evident  when,  in  Eq.  (76),  sin  <f>2  and  cos  $2  are 
replaced  by  the  expressions  given  in  Eqs.  (82)  and  (83): 


Ty  = 


(101) 
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l'z  =  -T||  sitl  *2  =~ 


sin 


21 


"*  rii  rn  rp  rp 

ui  =11  1  +  U  I  +  It  I 

o  x  x  y  y  z  z 


and  the  “penetration  depth,”  dp>  is  (Ref.  44) 


(102) 

(103) 

(104) 


dP  = 


to  n. 


sin2 


*  1 


n  * 


-  1 


‘21 


=  ^  2  ?r  nj  ^sin2  <^j  -  n2  t  ^  J 


-l 


(105) 


Equation  (99)  is  that  of  a  periodic  disturbance  that  is  propagating  in  the  rarer  medium  along 
the  interface  in  the  plane  of  incidence  (the  x-direction),  and  whose  amplitude  experiences  an 
exponential  decay  along  the  normal  away  from  the  interface  (the  z-direction). 


The  penetration  depth  is,  then,  infinite  at  the  critical  angle  tfI2c  =  sin  “  ln2i  -  The  smallest 
value  of  dp  is  at  grazing  incidence  (</>i  =  tt/2),  and  we  can  estimate  its  value  by  examining 
Eq.  (105)  at  tj>i  =  tt/2  and  for  very  small  n2i.  We  have 

£im  d  (0j  =  „/2)  =  (2  i/  n ,  )_1  =  - - -  (106) 

n2l-*a  2  ff  nj 

So,  for  two  transparent  media,  the  depth  of  penetration  can  be  anywhere  from  a  fraction, 
(2tt)  _  1 ,  of  the  wavelength  in  the  dense  medium  (X/ni)  to  infinite;  hence,  it  is  possible  to 
selectively  sample  a  thin  layer  in  multiple-layer  films.  Also  note  that  at  a  fixed  angle  of 
incidence,  the  penetration  depth  increases  with  wavelength. 

Correspondingly,  a  consideration  of  absorbing  media  requires  knowledge  of  the  relative 
magnitudes  of  the  transmitted  electric  fields  very  near  the  interface.  Expressions  are  derived 
herein  for  those  quantities. 


If  Eqs.  (80),  (81),  (86),  and  (87)  are  combined  with  Eqs.  (100)  through  (102),  we  may 
write 
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(107) 


(108) 


(109) 


(HO) 


(111) 


(112) 
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Also, 


The  factor  of  2  in  each  of  these  four  relations  indicates  that  the  evanescent  field  is  in  the 
form  of  a  standing  wave  (Ref.  44). 

3.1.2  Attenuated  Total  Reflection  for  a  Semi-Infinite  Transparent  Dense 
Medium  and  a  Semi-Infinite  Absorbing  Medium 

If  the  transparent  rarer  medium  of  the  last  section  is  replaced  by  a  semi-infinite 
absorbing  rarer  medium,  there  is  energy  coupling  (ATR)  by  the  incident  beam  with  the 
evanescent  field  in  the  rarer  medium  when  the  incident  angle  exceeds  the  critical  angle. 
Actually,  because  of  the  limited  extent  of  the  evanescent  wave  into  the  rarer  medium,  it  is 
usually  unnecessary  to  experimentally  approximate  a  semi-infinite  with  a  large  block  of  rarer 
material  volume.  A  thick  film  of  the  rarer  medium  will  suffice  if  its  thickness,  d,  is  much 
greater  than  the  penetration  depth,  dp. 

In  the  case  of  an  absorbing  rarer  medium,  with  complex  index  of  refraction 
n2  =  n2  +  ik2,  the  critical  angle,  (j>l2c,  is  not  given  by  Eq.  (72)  because  fi2  is  not  real. 
Furthermore,  the  amplitude  of  a  wave  propagating  within  an  absorbing  medium  depends  on 
the  distance  it  travels  within  the  medium. 

Internal  reflection  measurements  of  the  attentuated  total  reflection  processes  at  the 
interface  of  a  transparent  medium  (with  an  index  of  refraction  nt  =  n3)  and  an  absorbing 
medium  (with  a  complex  refractive  index  of  n2  =  n2  +  ik2)  can  be  interpreted  from  the 
complex  Fresnel  relations  [Eqs.  (58)  through  (61)]  when  modified  by  Snell’s  law  [Eq.  (70)] 
and  complex  arguments.  Thus,  at  the  interface  of  media  i  and  j, 

cos  =  yj l  -  sin2  4>.  =  yj  1  -  (njn.)2  sinS  (114) 
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the  complex  Fresnel-Snel!  relations  are 


rup 


n.  -\/  1  —  (n. /n.  )2  sin2  —  n.  cos 
t  \ _ i  /  i _ j _ J _ i_ 

n."^l  —  (n./ii.  )2  s'"2  +  "j  cos 


(115) 


2  n.  cos  <t>. 


and 


l‘jp 


ljS 


l01s 


a:\f  1  “  (n./n  sin^  4>.  -t-  n.  cos  <t>. 

1  V  J  i  j  i 

n.  cos  $.  —  n 1  ~  (iij/n.  )2  sin2 
n.  cos  4>.  -t-  1  ~  (n.^n.  )Z  sin2  4>. 

2  n.  cos  <5. 


n.  cos  <f> ^  +  n.  yj  1  —  (n.^n.  )"  sin^  0. 


(116) 


(117) 


(118) 


These  four  relations  apply  to  all  TIR  and  ATR  phenomena  at  the  interface  of  sufficiently 
thick  pieces  of  optically  dense  and  rare  materials,  whether  or  not  the  media  are  absorbing. 


The  inversion  of  these  formulas  to  obtain  the  optical  constants  usually  cannot  be  done 
exactly,  and  a  computer  algorithm  such  as  Marquardt’s  (Ref.  34)  must  therefore  be  used. 
Approximate  solutions  for  the  optical  properties  of  weakly  absorbing  media  are  useful  and 
easily  obtainable,  however,  using  the  concept  of  effective  thickness  (Ref.  44).  These  are 
found  in  the  next  section. 


3.1.3  Effective  Thicknesses  of  a  Semi-Infinite  Weakly  Absorbing  Rarer  Medium 

If  the  semi-infinite  incident  medium  is  transparent,  and  the  semi-infinite,  or  thick 
(d  >  >  dp),  rarer  medium  is  nearly  transparent,  then  to  a  good  approximation,  the  critical 
angle,  0i2c,  is  that  of  Eq.  (72).  The  effective  penetration  depth,  de,  is  defined  from  the 
relation 

where  a  =  4jrkv  is  the  Lambert  absorption  coefficient  in  the  Lambert-de  Beers  law  of 
intensity  transmission  in  the  z-direction: 

dl  =  -I  adz,  or  I  =  10  e‘az  (120) 


for  intensity,  I  (Ref.  42). 
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Harrick  (Ref.  47)  has  considered  the  low-absorption  rarer  medium  as  a  perturbation  of 
the  nonabsorbing  rarer  medium  situation;  he  gives  for  the  effective  penetration  thickness 


n  | T/E | 2  d 

d  =  21 

e 


2  cos 


A  /V 

n2 1  T/E 


4  it  v  n  t  ^ain2  —  r2  ^  cosdj 


(121) 


where  Eq.  (105)  is  substituted.  Using  either  Eq.  (Ill)  or  Eq.  (113)  in  Eq.  (121)  produces 


n^j  cos 


01  (2  3ii2  <£t  ~n|x) 


w»'nl(1  ~n2l)  (sin2  01  -n2l)A(si"2  <*1  -n21  cos2  01  ) 


(122) 


or 


n2  1  coa  0t 


nv  n](2  ~ali)  (si"20, 


(123) 


Obviously,  an  increase  in  the  effective  penetration  depths  can  result  from  decreasing  </>i 
and/or  v,  and  better  refractive  index  matching  (i.e.,  letting  n2i  approach  1).  Note  that  de|  is 
greater  than  dex  at  a  given  angle  of  incidence;  therefore,  the  parallel  component  alone  is 
frequently  used  in  experiments  on  a  thick  film. 

3.1.4  Effective  Thicknessess  of  a  Very  Thin  Weakly  Absorbing  Film 

In  this  section,  a  very  thin  film  is  considered  to  be  one  whose  thickness,  d,  is  very  much 
smaller  than  the  penetration  depth  computed  for  two  transparent  media  [Eq.  (105)], 
d  <  <  dpo.  Let  the  very  thin  film  be  isotropic,  planar,  and  of  complex  refractive  index,  n2. 
The  thin  film  is  sandwiched  between  two  semi-infinite  transparent  media:  the  incident 
medium  ni  and  the  medium  nj.  V 

At  the  interface  of  media  n;  and  nj,  the  fraction  of  the  incident  power  entering  the  second 
medium  nj  is  |  nj/nj  1 1  tjj  | 2  =  (n/nj)  |  tjj  | 2  for  transparent  or  nearly  transparent  media.  If  the 
power  is  assumed  to  be  constant  with  film  thickness,  then  Eqs.  (Ill)  and  (113)  must  be 
modified  to  read  (Ref.  42) 
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A  j 

Tj^  |  2  cos  <p  j 


(124) 


and 


IL 


2  cos 


4>1  [(*  +  n32)sinZ  jl  ~  nSl] 


where 


'll  {(!  -nL)(3in2^l  "n31  COs2^l)} 

"Jl  "  ni/"i 


For  the  very  thin  film,  the  effective  thickness  is  (Ref.  42) 

T/E  I2  d 


de  = 


21 


COS 


which  becomes,  upon  combining  with  either  Eq.  (124)  or  (125),  either 

4  J  Hj!  cos 


°1 


or 


tdn^cos^  f  (l  -  n32)sir*2  n3l] 

{  (  1  ”  n3l)  (si"2  *1  “  "31  cos2  ^l)}^ 


(125) 


(126) 


(127) 


(128) 


(129) 


Here,  de.  may  be  greater  or  smaller  than  dei. 


3.2  Experimental  Techniques 

The  two  major  divisions  of  internal  reflection  spectroscopy  are  single  reflection  methods, 
or  ATR  methods,  and  multiple  internal  reflection,  or  MIR  techniques.  For  both  kinds  of 
internal  reflection  spectroscopy,  the  proper  internal  reflection  element  (IRE)  is  very 
important. 
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3.2.1  Internal  Reflection  Elements 

In  Ref.  42,  Harrick  discusses  single  and  multiple  IRE’s  at  length.  He  notes  that  the  IRE 
material  must  have  a  high,  known  index  of  refraction,  n,  and  must  be  transparent  or  very 
nearly  so  (especially  for  MIR  measurements)  because  of  the  possible  loss  of  beam  power 
when  the  beam  is  within  the  IRE. 

Although  fixed-angle  IRE’s  (such  as  prisms)  are  usually  simpler  to  machine,  variable- 
angle  IRE's  (such  as  hemicylinders)  are  necessary  if  measurements  are  desired  at  several 
angles  (Ref.  40).  For  weakly  absorbing  samples,  the  number  of  multiple  absorptions  can  be 
increased  by  designs  involving  multiple  passes  of  the  beam  through  the  IRE. 

While  there  are  many  IRE  materials  (such  as  Ge,  KRS-5,  and  NaCl)  whose  optical 
properties  are  known  in  the  visible  and  IR  regions  at  300°K,  very  few  optical  constants  are 
known  for  these  or  other  suitable  materials  at  cryogenic  temperatures.  Ge  has  been 
measured  at  20°K  and  80°K  (Ref.  28),  but  other  cryogenic  measurements,  especially  on 
KRS-5,  should  be  undertaken. 


3.2.2  Single  Reflection  Methods 

If  the  samples  have  sufficient  absorption,  a  single  reflection  technique  is  often  simpler 
and  preferable  to  the  multiple  reflection  techniques.  Many  investigators  have  applied  single 
reflection  techniques  to  studies  of  thick  solids  and  liquids  throughout  the  IR  (Refs.  40,  42, 
46,  47,  48,  and  49). 


In  most  single  reflection  studies  the  incident  beam  is  unpolarized,  but  the  spectrometer 
optics  introduce  a  partial  polarization.  If  the  spectrometer  passes  a  ratio,  S(v),  of 
perpendicular,  s,  to  parallel,  p,  intensity  at  the  wavenumber,  v ,  then  the  measured  quantity, 
R,  is  the  ratio  of  the  beam  intensity,  I,  when  a  sample  is  in  contact  with  the  IRE,  to  the  beam 
intensity,  IOI  when  no  sample  is  on  the  IRE.  Let  nj,  n2,  and  nj  refer,  respectively,  to  the  IRE 
material,  the  sample,  and  the  ambient.  The  following  results,  then,  for  thick  samples 

"  -  |!I2|2/M!  (I30) 

where 


i  +■  s 


(131) 
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and  ryi  is  obtained  from  Eqs.  (1 15)  and  (116).  For  a  prism  instrument,  S  is  nearly  constant 
over  a  typical  absorption  band,  whereas  a  grating  instrument  may  have  a  wildly  varying  S 
over  the  same  spectral  region  (Ref.  46), 

For  thick  samples,  Eqs.  (115),  (116),  and  (131)  show  that  R  is  a  function  of  nine 
parameters:  the  wavenumber,  v;  the  angle  of  incidence;  4m;  the  instrument  polarization 
function,  S;  and  the  real  and  imaginary  parts  of  n1(  n2,  and  n3.  In  principle,  with  at  least  nine 
independent  measurements  of  R,  one  could  find  the  values  for  all  nine  parameters  at  the 
computed  wavenumber.  Usually,  k]  =  k3  =  0(i.e.,  media  1  and  3  are  not  absorbing),  and 
m,  n3,  v,  4m,  and  S  are  measured,  thereby  reducing  the  number  of  unknown  parameters  to 
the  two  optical  constants,  n,  and  kt>  of  the  thick  sample.  Inversion  of  Eq.  (130)  can  only  be 
done  iteratively  by  a  computer. 

With  a  single  hemicylinder  IRE,  Crawford  et  al.  (Ref.  46)  first  measured  S(v)  over  the 
fairly  small  wavenumber  region  of  interest  (typically  about  200  cm-1),  using  a  stack  of 
inclined  AgCl  plates.  Using  a  wide  range  of  incident  angles,  they  then  obtained  a  preliminary 
determination  of  the  optical  constants.  They  finally  used  a  program  which  computes  an 
optimum  set  of  angles  to  be  used  in  the  final  measurement  of  the  optical  constants. 

On  the  other  hand,  Irons  and  Thompson  (Ref.  49)  used  several  prism  IRE’s  made  of 
different  materials;  the  prisms  had  different  fixed  angles  of  incidence.  The  optimum  angles 
of  incidence  were  determined  empirically  for  each  prism,  and  several  values  of  R  were  found 
at  a  given  wavelength . 

3.2.3  Multiple  Internal  Reflection  (MIR)  Techniques 

MIR  techniques  are  preferable  to  single  reflection  methods  for  spectral  regions  in  which 
the  absorption  of  thick  materials  is  weak,  or  for  such  cases  when,  the  sample  for  such  cases 
being  a  very  thin  film,  the  absorption  of  energy  per  reflection  from  an  incident  beam  is  too 
small  to  measure  accurately. 

If  a  fraction,  R,  of  the  incident  intensity  is  reflected  per  reflection,  then  the  measured 
quantity  in  an  MIR  experiment  with  N  reflections  is  the  ratio  of  the  beam  intensity,  I,  when  a 
sample  is  in  contact  with  the  IRE,  to  the  beam  intensity,  Iot  when  no  sample  is  on  the  IRE. 
Thus, 


Inversion  of  Eq.  (132)  requires  a  computer  algorithm. 
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In  what  follows,  ni  refers  to  the  IRE  (dense)  material,  nj  to  the  sample  and  113  to  the 
ambient  index. 

3.2.3. 1  MIR  Methods  for  Thick  Samples 

For  samples  whose  minimum  thickness,  d,  is  much  larger  than  the  penetration  depth,  dp 
[Eq.  (105)],  R  is  the  same  as  given  in  Eq.  (130)  for  all  thick  samples.  If  the  radiation  is  plane- 
polarized  parallel  to  the  plane  of  incidence  (i.e.,  S  =  0),  then  R  =  |ri2pl2/|ri3p|2»  and  if 
polarized  perpendicularly  (i.e. ,  S  =  0),  then  R  =  |  r12p  | 2/ 1  rup  | 2.  The  measurements  I/l0  are 
functions  of  10  parameters:  v,  <j>\ ,  S,  N,  and  the  real  and  imaginary  parts  of  nj ,  n2,  and  113,  of 

which  all  but  n2  and  k2  of  the  sample  are  usually  measured  before  an  MIR  experiment. 

/ 

The  use  of  MIR  methods  for  solid  and  liquid  thick  samples  is  discussed  by  Rawlins  (Ref. 
48),  Medeck  (Ref.  50),  and  others.  (More  extensive  lists  of  studies  are  given  in  the 
bibliographies  of  Refs.  42  and  43.) 

3.2.3. 2  MIR  Methods  for  Thin  Films 


If  the  penetration  depth,  dp,  is  greater  or  comparable  to  the  film  thickness,  d,  then  the 
proper  value  for  R  is 


Rfilm 


(133) 


where 


12Sp 


123s 


123 


1  +  S 


(134) 


and  rmi  is  analogous  to  r<>]2i  >n  Eq.  (67).  For  polarized  radiation,  Rnim  is  found  from  the 
appropriate  component  values  of  rm  and  r^  inserted  into  Eq.  (133).  The  same  10 
parameters  appear  in  Eqs.  (130)  and  (133).  As  in  the  thick  sample  case,  though,  usually  all 
but  the  sample  optical  constants  m  and  ki  are  known.  Again,  a  computer  algorithm  such  as 
that  used  by  Marquardt  (Ref.  34)  is  required  to  solve  for  any  unknown  parameters. 


If  the  optical  constants  of  thin  films  are  to  be  determined  by  MIR,  preliminary 
calculations  may  be  made  using  the  effective  depths,  given  in  Section  3.1.4,  to  find  tentative 
values  of  the  critical  angle,  an  appropriate  IRE  material,  the  proper  geometry  of  the  IRE, 
and  other  experimental  parameters.  These  effective  depths  may  also  be  used  in  certain  cases 
to  obtain  good  approximate  values  of  the  optical  constants  of  the  sample  (Ref.  44).  The 
accurate  determination  of  the  optical  constants  may  require  several  experimental  runs, 
however,  as  well  as  an  analysis  of  the  data  using  Eq.  (133). 
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\  . 


Numerous  authors  have  used  MIR  spectra  of  thin  films  for  qualitative  analyses,  rather 
than  for  finding  the  optical  properties  of  the  sample.  (See  Refs.  41,  42,  44,  51,  and  52.) 

3.3  COMMERCIAL  INSTRUMENTATION  FOR  INTERNAL  REFLECTION 

SPECTROSCOPY  (IRS) 

Instrumentation  for  IRS  is  available  from: 

1.  Barnes  Engineering 
Stamford,  Connecticut 

2.  Beckman  Instruments,  Inc. 

Fullerton,  California 

3.  Perkin-Elmer  Corp. 

Norwalk,  Connecticut 

4.  Research  and  Industrial  Instruments  Co. 

London,  England 

5.  Foxboro  Analytical 
Wilks  Infrared  Center 
South  Norwalk,  Connecticut 

The  pioneering  manufacturers  were  Connecticut  Instruments,  which  built  the  first 
commercial  IRS  attachment  for  spectrometers  in  1965  (and  is  now  a  part  of  Barnes 
Engineering)  and  Wilks  Scientific  Corp.  (now  a  part  of  the  Foxboro  Company).  The  first 
IRS  units  were  accessories  for  existing  sing-  and  double-beam  spectrometers  including  single 
reflection  and  MIR  units.  There  are  now  IRS  units  available  for  microsamples  (with  a 
diamond  IRE),  for  variable-angle  MIR  spectroscopy,  and  for  skin  analyses,  and  narrow 
MIR  plates  made  of  KRS-5,  Ge,  Si  and  ZnSe.  There  are  now  specialized  IRS  spectrometers 
that  can  make  routine  analyses  rapidly  with  microcomputer  control. 

It  is  difficult  to  rate  the  commercial  IRS  products  available.  Most  IRS  units  but  not  all 
can  be  attached  to  the  popular  IR  spectrometers  including  Fourier  Transform  instruments. 
Some  units,  however,  may  be  more  easily  placed  in  position  than  others,  depending  on  the 
spectrometer.  As  a  general  rule,  the  units  designed  by  a  manufacturer  for  its  own 
spectrometers  are  more  easily  used  with  those  instruments.  The  authors  have  noted  that  in 
the  articles  they  have  read,  Wilks’  instrumentation  was  utilized  more  than  other  equipment. 
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3.4  SUMMARY  ON  INTERNAL  REFLECTION  SPECTROSCOPY  OF  FILMS 

In  the  1R,  visible,  and  near-UV  spectral  regions,  IRS  spectra  may  often  be  superior  to 
conventional  reflection  or  transmission  spectra,  especially  for  weakly  absorbing  samples  or 
very  thin  films.  IRS  techniques  may  be  used  in  conjunction  with  ellipsometric  measurements 
on  thin  films  because  the  IRS  measurements  are  sensitive  to  the  optical  constants  of  the 
material,  and  ellipsometer  data  depend  heavily  on  film  thickness.  Films  having  a  minimum 
thickness  much  greater  than  the  penetration  depth  of  the  evanescent  field  have  IRS  spectra 
which  are  free  of  multiple  reflection  effects,  even  for  nonplanar  films.  Accurate  reduction  of 
IRS  data  for  films  of  any  thickness  can  be  done  only  by  computer  programs. 

Finally,  a  wide  variety  of  stock  items  for  IRS  is  currently  available  for  single  reflection 
and  MIR  studies  of  liquids  and  solids. 

4.0  ANALYSES  OF  MOLECULAR  ABSORPTION  BANDS  BY 
THE  METHOD  OF  MOMENTS 

A  crucial  part  of  accurate  determination  of  mole  fractions  of  IR-active  species  in 
cryodeposits  is  the  calculation  of  band  strengths  of  those  species.  Discrepancies  noted  earlier 
(Refs.  1  and  2)  between  calculated  mole  fractions  and  those  obtained  by  chemical  analysis 
were  thought  to  be  partly  caused  by  experimental  uncertainties  in  the  computation  of  band 
strengths.  In  particular,  choices  of  integration  limits  and  background  absorption  (baseline) 
were  made  on  an  intuitive  or  empirical  rather  than  quantitative  basis.  It  was  recommended 
at  that  time  that  a  “moment”  method  be  tried  and  evaluated  in  future  band  strength 
calculations,  as  a  possible  means  of  reducing  experimental  uncertainty. 

In  a  series  of  papers  (Refs.  54  through  57),  Gordon  analyzed  vibration  bands 
theoretically  in  terms  of  their  moments.  By  using  data  obtained  from  various  experimental 
studies,  he  tested  the  theoretical  analysis  for  consistency  and,  in  the  process,  provided  a 
method  to  determine  integration  limits  and  baselines  quantitatively.  (For  related  work  on 
band  analyses,  see  Refs.  58  through  66.) 

4.1  THE  FIRST  AND  SECOND  MOMENTS 

Much  of  the  present  work  was  devoted  to  moment  analysis.  The  scheme  adopted  follows 
that  suggested  by  Gordon’s  work  and  implemented  by  Jiang  et  al.  (Ref.  67).  The  moment 
method  rests  on  the  assumption  that  the  first  and  second  moments  of  an  absorption  band 
are  independent  of  intermolecular  forces.  Thus,  the  calculated  second  moment  must  match 
the  “classical”  second  moment  plus  quantum-mechanical  correction  terms  that  become 
significant  at  low  temperatures. 


54 


A  E  DC-TR  -80-51 


In  general,  the  nth  moment  is  defined  (Ref.  67)  by 

M(n)  =  /(a»  -  )  p  (<y)  dw  jj p  (<w)  du>  (135) 

1 

where  w5  is  the  “shifted  band  origin,”  p(u),  is  an  intensity  distribution  function,  and  oi  = 
2kcv.  The  shifted  origin,  os,  is  chosen  so  that  M(l)  agrees  with  the  classical  first  moment  Mi: 


MU)  = 


f  (to  —  G>a  )  P  (dj)  dcu 

/  p  (cj)  dw 


(136) 


4.2  COMPUTER  PROGRAM  MOMENT 


A  moment  analysis  program  has  been  written  using  moments  converted  to  v  units  with 
the  distribution  function  p(w)  =  p(2ircp)  =  k(p).  The  computer  program  MOMENT 
contains  five  options  for  k(p).  AEDC  k-values  were  used  in  addition  to  those  k-values  fitted 
by  nonlinear  least  squares  to  the  following  four  models: 

kc(v)  -  - 2 - 

1  +  c2  U--b)2 

kg  (.v)  ~  F  ■  exp  £-d2  (v  -  b)2J 

kpM  =  kc  w-kgw 

k,w  =  kcMtljM  (137) 

The  models  kc  and  kg  are  the  familiar  Cauchy  (Lorentz)  and  Gauss  profiles;  they  provide  one 
possible  method  of  treating  blended  absorption  bands  (bands  such  that  the  band  of  interest 
is  blended  with  a  lattice  vibration).  Jones  et  al.  (Refs.  69  through  71)  have  used  nonlinear 
least  squares  to  fit  absorption  spectra  and,  in  some  instances,  have  found  the  product  or  sum 
model  to  give  a  better  fit  than  the  pure  Cauchy  or  Gauss  profiles. 


Program  MOMENT  was  designed  to  first  calculate  the  shifted  band  origin,  vs,  using  Eq. 
(136),  and  to  then  force  the  calculated  second  moment,  M(2),  to  match  the  classical  second 
moment,  M2i  plus  quantum  corrections.  All  calculated  moments,  M(n),  n  >  2  are  computed 
with  respect  to  ps,  and  the  equation  used  is  that  derived  by  Gordon  (Ref.  58): 


M(2)  =  M,  +  -HilL  .  (hi  V  _  /i-V  +  17MW  /he  Y  - 
12  \kT/  120  \kT/  20,100  \  kTV  ' 


038) 


The  right  side  of  Eq.  (138)  is  an  infinite  alternating  series,  convergent  whenever  |i»]  < 
irkT/hc  *  43.67  and  asymptotic  outside  that  interval.  An  iterative  scheme  was  used  in  an 
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attempt  to  force  the  quality  indicated  in  Eq.  (138),  by  using  corrections  through  the  M(8) 
term.  The  (horizontal)  baseline  height  and  limits  of  integration  were  adjusted 
(simultaneously,  in  the  case  of  the  model  options)  for  SO  iterations  or  until  agreement  of  the 
order  10“ 5  was  reached.  The  limits  of  integration  were  chosen  to  be  those  ^-values  for  which 
the  baseline  intersected  the  k-profile.  When  the  AEDC  k-values  were  used,  Simpson’s  rule 
(modified  for  unequal  subdivisions)  was  used  to  calculate  the  moments;  when  model  options 
were  used,  the  integrations  were  performed  either  analytically  or  by  using  10-point  Gaussian 
quadrature  over  each  of  10  equal  subintervals  of  the  integration  range. 


4.3  ANALYSES  WITH  PROGRAM  MOMENT 

The  MOMENT  program  was  first  used  to  analyze  p2  of  H20,  l2CO(l  —0),  and  vj  of 
”C02;  AEDC  data  from  20°K  cryofilms  were  used.  The  iterative  scheme  failed  to  converge 
in  all  cases  because  the  calculated  second  moment,  M(2),  was  so  small  compared  to  the 
classical  value,  M2.  Either  truncation  of  the  wings  or  insufficient  resolution  of  data  in  the 
wings  could  cause  these  weak  second  moments. 

All  of  the  models  were  then  fitted  by  nonlinear  least  squares  to  the  bands  mentioned  in 
the  previous  paragraph;  only  one  positive  result  was  obtained,  that  with  the  Cauchy  fit  to  *3 
of  !3C02.  Again,  in  all  but  this  one  case,  the  calculated  second  moment  was  well  below  the 
classical  value,  or,  when  the  baseline  and  limits  were  adjusted  to  achieve  a  closer  agreement 
between  M(2)  and  M2,  the  correction  terms  became  extremely  large,  and  the  correction  series 
became  divergent.  Since  accurate  representation  of  absorption  bands  in  the  wings  is  essential 
for  successful  moment  analysis,  the  models  might  require  modification  or  other  models 
might  have  to  be  found  that  better  approximate  bands  in  the  wings. 

Incidentally,  the  computer  results,  even  in  the  failed  cases,  offered  empirical  evidence  of  1 
some  general  results  on  moment  analysis  with  symmetric  models,  some  of  which  may  be 
used  in  practice  to  determine  whether  a  given  model  adequately  represents  k(p).  Some  of 
these  results  were  verified  mathematically  and  are  discussed  later  in  the  report. 

By  using  the  Cauchy  representations  of  v±  of  l3C02  in  pure  C02  and  in  the  binary 
mixture  N2/C02,  new  values  were  obtained  for  the  band  strengths  in  both  the  standard 
cryofilm  (pure  C02)  and  in  the  mixture.  The  moment  analysis  produced  “corrected” 
baselines  and  limits  of  integration  which  gave  excellent  agreement  between  M(2)  and  the 
quantum-corrected  M2.  The  difference  M(2)  -  M2  was  less  than  10s  in  both  cases.  The  band 
strength  ratio,  R,  and  mole  fraction,  f,  were  calculated  and  compared  with  previously 
reported  values  (Ref.  1).  The  previous  value  obtained,  f  =  0.296,  is  within  17  percent  of  the 
mole  fraction  determined  by  chemical  analysis.  The  present  value  obtained,  f  =  0.289 
represents  a  slight  improvement,  falling  within  14.3  percent  of  the  chemically  determined 
value  of  0.253. 
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4.4  MATHEMATICAL  BAND  MODELS 

To  help  assess  the  appropriateness  of  the  various  mathematical  models,  it  is  useful  to 
examine  some  of  their  results.  It  is  assumed  that  each  distribution,  k(v),  is  symmetric  about 
the  band  center,  b,  and  that  the  interval  of  integration  is  (b  -  5,  b  +  5),  5  >  0.  (Detailed 
calculations  are  given  in  the  Appendix.) 

Result  1.  The  shifted  band  origin  is  independent  of  the  baseline  and  integration  limits.  In 
fact,  b  -  i>s  =  Mj,  the  classical  first  moment. 

Result  2.  For  the  Cauchy  profile, 

Sim  M  (n)  =  do  for  n  2 

8  -*  DO 

as  Gordon  and  others  have  pointed  out. 

Remark.  Curiously,  it  is  this  property  of  kc(p)  that  makes  it  possible  for  our  moment 
analysis  to  succeed  in  those  cases  when  M{2)  -  M2  is  sufficiently  small  and  positive. 

Result  3.  For  the  Gaussian  profile,  kg(i>),  all  moments  are  finite  as  6  —  <*.  In  particular, 

fim  M (2)  =  (l/2d2)  +  M2  (139) 

Here  d  =  2(ln2)l/z/W,  where  W  is  the  band  width  at  half  height. 

Remarks.  Equation  (139)  can  be  used  in  practice  to  decide  whether  the  Gaussian  profile 
is  appropriate  for  a  given  band.  For  example,  the  formula  in  Eq.  (139)  yields  an  upper 
bound  of  5.0052  for  M(2),  using  the  vj  band  of  CO2,  (In  this  analysis,  parameters  such  as  d 
are  converged  least-squares  values  for  the  Gaussian  fit  to  the  AEDC  data.)  The  classical 
second  moment  for  C02  is  M2  =  21 .862858.  Thus,  the  MOMENT  analysis  method  will  not 
succeed  for  this  Gaussian  profile,  regardless  of  how  the  baseline  and  limits  are  adjusted, 
since  M(2)  has  a  fixed  upper  bound  far  below  M2. 

Result  4.  The  higher  moments  are  calculated  about  the  shifted  band  origin,  with  baseline 
1  =  k(b  -  6)  =  k(b  +  6): 

M(n)  =  /(^-vs)n  [k(v) (l4°) 
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These  moments  can  be  rewritten  in  terms  of  moments  about  the  axis  of  symmetry,  v  =  b, 
where  b  is  the  band  center.  The  even  moments,  in  particular,  are  given  as: 


where 


M(2n)  = 


•  Mb <2n - 2P  - T j/u - 2*s) 


T  =  {1/ 2n  +  1)  [(5  +  M,)2n+  1  +  (S-Mj)2"" 


Mb(n)  =  /(v-b)n  k  (*)  dv 

and 

j  =  y  k  (^)  <1^  =  Mb  (o) 

Remark.  Writing  the  moments  in  this  manner  allows  the  integration  to  be  done 
analytically  for  the  models  in  Eq.  (137).  We  rewrote  the  moment  program  for  the  Cauchy 
model  on  a  Radio  Shack®  TRS-80  microcomputer,  doing  all  integrations  analytically,  and 
obtained  results  very  similar  to  those  of  the  original  MOMENT  program.  The  rewritten 
program  should  and,  we  believe,  did  give  slightly  better  values  for  the  moments  and  band 
integrals. 


4.5  COMPARISON  WITH  OTHER  MOMENT  ANALYSES 


All  of  the  above  analysis  dealt  with  20DK  cryofilms.  In  addition  to  these  attempts,  we 
tried  to  duplicate  the  results  of  Jiang  et  a).  (Ref.  69)  using  their  band  parameters  and  a 
Cauchy  model  for  the  band  12CO(l  —0)  at  30°K.  This  attempt  failed,  again  because  of  a  very 
weak  calculated  second  moment  compared  to  the  classical  value.  Thus,  it  may  be  instructive 
to  consider  some  of  the  differences  between  the  present  analysis  and  that  of  Jiang  et  al.  First, 
they  used  the  distribution  function  p(w)  =  fti(Io/I)  from  an  absorbance  measurement,  while 
the  present  research  involved  k(p)  fitted  by  a  Cauchy  model. 

Jiang  et  al.  integrated  geometrically  (band  integrals,  for  example,  were  calculated  either 
from  data  plots  or  directly  from  the  spectrum  using  a  polar  planimeter),  getting  best  results 
for  a  baseline  that  matched  the  intensity  in  the  wings  and  followed  "the  estimated  fringe 
pattern  for  a  hypothetical  thin  nonabsorbing  CO  film.”  The  model  used  in  the  present  study 
for  k(j»)  is  smooth  in  the  wings,  with  "channel  spectra”  eliminated.  Also,  horizontal 
baselines  were  tried  and  symmetric  limits  chosen.  The  only  limits  given  by  Jiang,  et  al.  were 
in  a  footnote;  they  were  those  limits  applied  to  an  Ar/CO  mixture  at  30°K.  These  limits  were 
not  quite  symmetric  about  the  band  origin,  thus  giving  a  slightly  different  shifted  band 
origin  from  the  present  results.  Incidentally,  better  results  were  achieved  in  the  present  study 
by  choosing  skewed  limits. 
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4.6  SUMMARY  ON  MOMENTS 

The  present  modest  success  with  moment  analysis  proves  that  the  technique  deserves 
further  examination  and  refinement  and  that  the  questions  implied  by  the  above  discussion 
warrant  investigation.  In  particular,  questions  remain  to  be  answered  concerning  which 
models  are  best  to  use  for  the  intensity  distribution  function  and  how  the  limits  of 
integration  should  be  determined  as  the  baseline  is  adjusted. 

The  independence  of  the  second  moment,  M(2)t  from  intermolecular  interactions  (Refs. 
55,  56,  and  57)  especially  warrants  scrutiny.  It  is  hoped  that  a  method  can  be  developed 
centering  in  this  characteristic  that  would  allow  spectroscopic  determinations  of  species 
concentration  in  a  cryofilm  regardless  of  which  species  are  in  the  film. 

5.0  SUMMARY 

We  have  investigated  the  applicability  of  the  ellipsometric  and  internal  reflection 
spectroscopy  (IRS)  methods  of  reflection  spectroscopy  to  measurements  of  the  optical 
properties  of  thin  films  and  substrates.  Ellipsometry  is  very  sensitive  to  the  condition  of 
sample  surfaces  and  film  structure  and  can  accurately  measure  optical  constants  of  films  and 
their  thicknesses,  even  in  the  case  of  multiple  films.  The  technique  can  be  used  to  make 
optical  measurements  of  substrates,  but  these  data  must  be  treated  with  reserve  because  of 
the  possibility  of  spurious  effects  caused  by  molecular  absorption.  Ellipsometric 
measurements  in  the  IR  beyond  3  pm  require  special  methods  and  instrumentation  and  are 
not  as  accurate  as  observations  in  the  visible  range. 

i 

On  the  contrary,  IRS  data  can  easily  be  obtained  throughout  the  IR  spectral  range  as  well 
as  the  visible.  The  penetration  of  the  evanescent  wave  makes  IRS  more  a  probe  of  the  bulk 
properties  of  thick  samples  rather  than  of  their  surface  characteristics.  By  prudent  choices  of 
the  beam  incident  angle  and,  thus,  of  the  penetration  depth,  one  may  obtain  spectra  of 
strong  or  weak  absorption  features  that  have  higher  resolution  than  do  transmission  spectra 
or  other  reflection  data.  Accurate  IRS  data  of  thinner  films  are  also  made  possible  if 
multiple  reflection  effects  with  the  film  are  accounted  for.  One  may  also  capitalize  on  the 
lower  power  loss  per  reflection  and  use  multiple  internal  reflection  (MIR)  schemes  for 
weakly  absorbing  (thick  or  thin)  samples. 

To  extract  accurate  optical  constants  and  film  thicknesses  from  either  of  the  above 
reflectance  methods,  researchers  must  perform  a  computer  inversion  of  the  appropriate 
equations.  Because  of  the  numerous  methods  of  ellipsometric  or  IRS  measurements,  local 
algorithms  are  necessary  for  handling  the  experimental  configuration. 


59 


AEDC-TR -80-51 


The  analysis  of  molecular  absorption  band  moments  can  lead  to  a  quantitative 
evaluation  of  the  intermolecular  interactions  of  molecules  in  condensed  phases.  These 
interactions  cause  spectral  distortions  that  make  molecular  species  concentration 
determinations  difficult.  A  preliminary  investigation  of  moment  analysis  shows  that 
symmetric  band  shapes  must  meet  some  simple  criteria  in  order  to  find  an  accurate  value  of 
the  band  strength.  Also,  the  importance  of  accurate  profiles  of  the  wings  of  an  absorption 
band  in  moment  analysis  must  be  emphasized.  It  is,  perhaps,  because  of  this  fact  that  higher- 
resolution  data  are  necessary  for  meaningful  moment  analysis.  It  is  also  hoped  that  the  lack 
of  dependence  of  the  second  moment  upon  the  intermolecular  interactions  can  lead  to  a  way 
of  finding  molecular  concentrations  in  mixtures  varying  widely  in  their  constituents. 

Other  areas  of  consideration  for  future  research  are: 

1 .  analysis  of  the  physics  of  quartz  crystal  microbalances  (QCM), 

2.  computer  programming  for  the  analysis  of  ellipsometric  and  IRS  data, 

3.  investigation  of  the  effect  of  multiple  reflections  in  a  dielectric  substrate  during 
ellipsometric  measurements, 

4.  continuation  of  moment  analysis  and  the  determination  of  molecular  species 
concentrations  in  cryofilms, 

5.  completion  of  the  computer  program  (now  in  its  final  debugging  stage)  for 
Ktamers-Kronig  (KK)  analysis  of  the  transmission  of  a  beam  through  a  single 
thickness  thin  film,  and 

6.  completion  of  the  computer  program  for  KK  analysis  of  the  reflection  of  a  beam 
of  a  single  thickness  of  a  thin  film. 
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APPENDIX 

Here  are  established  the  report’s  results  concerning  moment  analysis  with  symmetric 
models. 

Result  1.  The  shifted  band  origin,  v%,  is  calculated  as  follows: 

/U  —  vs)  k(v)  dv/  /k  U)  di/  =>  Mj 

Thus, 

Mj  =  ^  j'(i/  -b)  k  (id  Av  +  /(b  —  v^ )  k  (j/)  drj  j /  k  U)  Av 


b  — 


since  j(r  -b)  k(r)  dr  =  0.  This  follows  from  the  facts  that  k(r)  is  symmetric  about  v  -  b  and 
that  the  integration  interval  is  (b  -  6,  b  +  5)  for  6  >  0.  The  same  result  follows  if  k(r)  is 
replaced  by  k{r)  -  1,  since  this  latter  distribution  has  the  same  symmetry  properties  as  k(r), 
and  the  new  limits  of  integration  will  remain  symmetric  about  b. 


Result  2.  It  is  easy  to  see  that  the  moments  M(n),  n  a  2,  for  the  Cauchy  model  become 
infinite  as  the  interval  of  integration  becomes  infinite.  In  general,  for  any  symmetric  model 
k(r), 

MU)  =  /(v-rs)n[kU)-f]dv//[kU)-f]dl, 


If  J  =  J  k(r)dr,  then 


_  {/[fr-^+th-v.)]  "  Hrf  dr-[^  +  } 


M  (n)  = 

J  -  228 

Binomial  expansion  in  the  fust  term  of  the  numerator  yields 


M  (n)  = 


-V  v  k  U)  d^  -T  /{J -228) 


where 

T 
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Because  of  symmetry,  the  odd  moments  about  b  are  zero,  so 


P 

— 

M  (n)  - 

i  =  o  V 1 

.  ^  M2!  Mfc  (n  2j)  T 

r(J  -  2 C8) 


where  p  is  the  greatest  integer  in  n/2  and 

\  M  =  /(v -b)n  It  {v)  Av  ,  =  b-t/s 

For  the  Cauchy  model, 

Mb(j)  =  *  f  l,-b,j/[l-c*(,-b)*]d. 


/•cS 

c  8 


j  +  1  J-r.8  1  +  u2 


du 


For  j  >  2,  and  even, 


and  is  divergent  as  8  —  ». 


-Vi)  -  fS  du 

c‘  +  1  J0  ]  +  u2 


Now, 


Resuit  3.  We  derive  Eq.  (139).  From  the  general  calculation  above  in  Result  2, 
M  (2)  =  |  Mfc  (2)  -  Mj  J  -  (C/3)  [(Mt  +  8)3  +  (Mj  -S)3]|^/(J  -  2 18) 

f 


(Mj  -  5)3  +  (Mj -S)3  =  25(3MJt82) 


Therefore, 


M 


(2)  =  [Mb  (2)  +  M2  J  -  (25 C/3)  (3M2  +  S2)]  fu  -  2t8) 


(2)  -  2  £  53/  3 
J  -  2P5 


+  M; 


For  the  Gaussian  model. 


Mb  (2) 


=  f^8  (v-  b)2  -f  .e  ‘d2 


(141) 


68 


AEDC-TR-80-51 


If  u  =  d(r-b),  then 


Thus, 


/dg  (u2/d2)'Ne'U2du/<! 

•  l/d3  (-i8t  d‘  S\  J‘ 


•dS  , 
c-u  du 


!±  -d 2s2  J 
d2  2d2 


M  (2)  = 


j/2d2  -  (fS/d2  )  e‘d  8  ~  2£53/a  +  m2 


J  -  2 £5 


2  _  2 

J  -  2fge-d  S  -  (4/3)gdz33  +  2 

2d2  (J  -  2fS)  1 

Now,  fe-d2s2  =  t,  so  f—  0  exponentially  as  6—00,  and 


£im  M  (2)  =  -i-  +  M? 

5-oo  2d2 

Result  4.  This  result  is  established  above  in  deriving  Result  2.  From  Eq.  (141),  we  need 
only  replace  n  with  2n. 
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